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Two-Body Dirac equations of constraint dynamics provide a covariant framework to investigate 
the problem of highly relativistic quarks in meson bound states. This formalism eliminates automat- 
ically the problems of relative time and energy, leading to a covariant three dimensional formalism 
with the same number of degrees of freedom as appears in the corresponding nonrelativistic prob- 
lem. It provides bound state wave equations with the simplicity of the nonrelativistic Schrodinger 
equation. Here we begin important tests of the relativistic sixteen component wave function solu- 
tions obtained in a recent work on meson spectroscopy, extending a method developed previously for 
positronium decay into two photons. Preliminary to this we examine the positronium decay in the 
'^Po,2 states as well as the ^Sq. The two-gamma quarkonium decays that we investigate are for the 
'7c) 'yci XcOi Xc2) ■"""i ■''"2) o^ajand /2 mesons. Our results for the four charmonium states compare well 
with those from other quark models and show the particular importance of including all components 
of the wave function as well as strong and CM energy dependent potential effects on the norm and 
amplitude. The results for the tt", although off the experimental rate by 15%, is much closer than 
the usual expectations from a potential model. We conclude that the Two-Body Dirac equations 
lead to wave functions which provide good descriptions of the two-gamma decay amplitude and can 
be used with some confidence for other purposes. 

PACS numbers: 12.39.Ki,03.65.Pm,12.39.Pn 

I. INTRODUCTION 

Relativistic treatments of the two-body problem arise in many problems in particle and nuclear physics. Relativistic 
effects are important for composite systems with light quarks, in systems with large coupling strength, and in reactions 
of these composite objects. In recent years, there is much interest in the dissociation and the recombinations of the 
J/ 4' particle in hadron matter or in the quark-gluon plasma [1]. Reactions of the form 

J/^ + TT^D + D*, (1.1) 

provide useful information on the suppression or the enhancement of J/i/' in high-energy heavy-ion collisions and are 
relevant to the use of heavy quarkonium as a diagnostic tool for the quark-gluon plasma [2] . 

Previously, Wong, Barnes, and Swanson studied the above reactions using a non-relativistic model of the reacting 
composite objects including pions [1, 3]. While the results have been calibrated with the tth scattering phase shifts for 
the 1 = 2 S-wave channel, the use of the non-rclativistic formalism for pions with light constituents may be subject 
to question. One should examine the reaction process using a well tested relativistic formalism. The Two-Body Dirac 
equations (TBDE) of constraint dynamics has had successful applications to relativistic two body bound states in 
QED [4, 5], QCD [6, 7], and two body nucleon-nucleon scattering [9, 10]. But its relativistic extension [11] of the 
nonrelativistic four-body scattering formalism of Barnes and Swanson [12, 13] involves untested assumptions beyond 
the standard constraint formalism. The reaction process is sensitive to the spatial distribution of the reacting objects. 
It is thus important to have a sensitive test of the wave functions obtained in [7]. 

We perform this test in this paper by examining the application of the relativistic constraint formalism in the 
description of decays of mesons into two photons. In the next section we present a brief review of the constraint 
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formalism as it applies to quark- antiquark bound states. Part of the purpose of this review section is to outline some 
of the numerous tests made so far on the formalism. We give the Pauli forms of the Two-Body Dirac equations of 
constraint dynamics that wc used in [7] to describe the entire meson spectrum (exceptions being light quark isoscalars 
such as the 77, w,?/ and their orbital and radial excitations). We review those aspects of the formalism which give 
one confidence in the accurate accounts for all bound states from the excited states of bottomonium to the pion. We 
also list those aspects particularly related to the pion and its Goldstone boson behavior. The constraint formalism, 
unlike most of the other ones that purport to account for the entire meson spectrum, accounts well in standard 
perturbative approaches as well as in nonperturbative and numerical approaches for the QED bound state spectrum 
[5]. We emphasize the importance of this by showing the correlation between that agreement for the singlet and 
triplet positronium systems and the tt and p states, including the Goldstone behavior of the former. 

This connection to QED brings us in Section 3 to our treatment of the 2 7 decay of singlet positronium. In 
constraint dynamics, the Two Body Dirac equations lead to an analytic solution of the singlet states of positronium. 
For the singlet ground state the wave function is mildly singular. Standard formalisms [14] will fail with wave 
functions that are singular at the origin, 

r(e+e- ^ 27) = atotVe+ \iiOf . (1.2) 



Independent treatments by Crater [15] and Ackleh and Barnes [16] develop related (but distinctly different) ap- 
proaches for folding in the effect of the Yukawa fermion exchange mass, giving a smearing of the singularity over the 
corresponding Compton wave length. We give a brief review of the first of these approaches and how we extend it 
to include the effects of the full sixteen component Two-Body Dirac wave function. This extension does not have 
any significant effect on the ^S'o positronium decay rate. However, the effects on the decays of the more relativistic 
quark- antiquark systems is significant. 

We include in Section 3 technical aspects in which wc establish in the context of a 4 x 4 matrix wave function, more 
natural for use in the decay formalism of a particle-antiparticle system than the 16 component form, the relation 
between the sector of the full wave function used in the Pauli form of the bound state equations and the remaining 
sectors necessary for a complete description of the decay. We review our ^6*0 positronium decay results as well as 
those of our constraint approach for ^Pp and positronium decay. Finally wc present the results for the decay 
rates of the Jycj '7ci XcOi Xc2i """i ""aj ^27 ./2 mesons. We conclude in Section 4 with a discussion of our results and a 
comparison with other approaches. 



II. CONSTRAINT DYNAMICS AND MESON BOUND STATES 

A. Constraint Dynamics for Two Classical Spinless Particles 

Here we give a brief review of the highlights of the constraint approach serving also to introduce notations. Although 
Sazdjian has shown that the bound state equations of constraint dynamics are to be viewed as "quantum mechanical 
transforms" of the Bethe-Salpeter equation [17]- [19] the constraint approach to the two body problem has its origins 
in classical relativistic physics [20]- [25]. Our review here is base on [8] and [22]. Two free spinless particles are 
described by the mass shell constraints 

nl =pl+mi « 0,7i^ ^pj+mj « 0. (2.1) 

We introduce Poincare' invariant world scalar interactions (to display most simply the basic ideas) by 

mi mi + Si{x,pi,p2) = Mi{x,pi,p2), 

TO2 ^ 1712 + S2{x,Pi,P2) = M2ix,Pi,P2), 

X ~ xi ^ X2- (2.2) 
Kinematical constraints then become dynamical mass shell constraints: 

H° = pi + ml ^ pf + Mf = = pi + mi + <^>,{x,pi,p2). (2.3) 
Each constraint must be conserved, implying that the two constraints must be compatible 
« {^1,^2} 

d d 

= -(Pi +P2) • ^(^'2 + $1) - (pi -P2) • ^($2 - *l) + {'I'l,$2}. (2.4) 
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Its simplest solution is 

$1 = $2 = $(a;±,pi,P2) = (2.5) 
and requires abandoning x = xi ~ X2 in favor of 

<2± - {rr + P''Pn{^i-^2)u^rfl^{x,^X2)., 

^ -^=^ ; P^=pt + P^2 ; ^12^ - P = 0- (2.6) 

Thus we have a "third law" condition (2.5) of action and reaction plus a restriction on how the quasipotential 
may depend on relative separation. The invariant r defined below is the interparticle separation in the CM frame 
P=(1,0) 

= \/r2 in CM frame P = (1, 0) , (2.7) 

since <i — ^2 = in that frame. Relative time is thus controlled in a covariant way. Assume the two invariants 
Mi, i = 1,2 are simply functions of r and the CM energy 



w = V^P2. (2.8) 

The invariant potentials Mi are not independent. The third law condition implies they are related by 

M^-M^ = ml-ml. (2.9) 

Hence there is only one independent invariant function controlling the scalar interaction which we designate by 

Sir), (2.10) 

the underlying scalar interaction. Alternatively, the third law allows us to recast the mass potentials into the hyperbolic 
function solutions depending on a single invariant function L, 

Ml = miCOshL{S{r))+ m2Smh L{S{r)), 

M2 = m2 cosh L{S (r)) +mismhL{S (r)). (2.11) 

Subtracting the constraints gives us a complimentary covariant restriction (to Eq. (2.7)) on the relative energy 

ni-n2^pl + Mf -pI- Ml ^pI + ml -pI- m\ = 2P ■ p wO, (2.12) 

with relative momentum 

P^ = , ei+e2 = w , ei - £2 = 

w w 

Ei = CM energy of particle i. (2.13) 
The relative momentum is canonically conjugate to x^, 

= (2.14) 
The other combination of our constraints is the primary dynamical equation 

n ^ =pI+ - b\w) « 0, (2.15) 
w 



and incorporates exact two-body kinematics with 

h\w) = i^'-'^^'irnl + ml) + [ml-ml)-) ^ _ ^^.16) 



and 

mim2 [w^ — m\ — m|) 



(2.17) 



defined as the mass and energy of the fictitious particle of relative motion. Under quantization all of the constraints 
become equations the wave functions must satisfy. 
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B. Two-Body Dirac Equations 

The constraint formalism embodies spin in a system of two coupled, compatible Dirac equations on a single wave- 
function. For particles interacting through world vector and scalar interactions the TBDE take this general minimal- 
coupling form 

Sii^ = 75i(7i ■ (Pi - ^i) +mi+ Si)i' = 0, 

S2i^ = 752(72 ■ (P2-i2)+ ma + 52)^' = 0. (2.18) 
The wave function has sixteen components 

V'= [^l,V2,^3:V'4], (2.19) 

in which each ■0^ is a four component Pauli spinor for two spin-one-half particles. The two equations are compatible: 

[5i,52]V = 0. (2.20) 

This is a result of the presence of spin supersymmetries [8], [26], in addition to the relativistic third law, and covariant 
restrictions on the relative time and energy appearing in the spinless case. There is automatic incorporation of correct 
spin-dependent recoil terms [27], 

if = if (A(r),pj., P, w, 7i, 72), = S,{S{r), Air),p^ , P, w, 71, 72)- (2-21) 

This two-body formalism has many advantages over the traditional Bethe-Salpeter equation and its numerous three 
dimensional truncations. One is its simplicity. A Pauli reduction and scale transformation brings our equations to 
this covariant Schrodinger-like form 

+ $^(ai, (72,p±, ^(r), S{r))^ i' = b\w)^. (2.22) 



1. Schro dinger- Like Form of the Two-Body Dirac Equations 

From separate classical [28] or quantum field theories [19] one can show that the quasipotential in the combination 
'^tu — (w) depends on the difference of squares of the invariant mass and energy potentials {Mi and Ei respectively) 

Mf = m\ + 2m.^S -f S^; = ef - 2e^A + , (2.23) 

with A playing the same role for vector interactions that S does for scalar ones. "Squaring" the TBDE (2.18) yields 
a Schrodinger-like equation [5] for the upper- upper V-'i component 



{p2 + 2m.^S + 52 + 2eu,A - A^ 

+ $_Dif • p + <^D' + '^SOlL ■ <Ji -\- ^S02L ■ Gi -I- <i>SSCri • CT2 -|- $T'S't}V'i 

= (2.24) 
coupled to a Schrodinger-like wave equation for the lower-lower component "04 [29] 

{p2 + 2m,„5 5*2 + 2e,,A - A^ 

-|-$_Dif • p + + "SsOli • fi + $502-^ • 02 + 'fsSO'l • CT2 + '^tSt^'^/^ 

+{^sscri ■ 02 + $t'5't}V'i 
= \?(w)^^. (2.25) 

These equations can be solved nonperturbatively for QED [S = 0) or quark model calculations since everyone of 
the quasipotentials terms (including the Darwin pieces $13) is quantum mechanically well defined (less singular 
than -l/4r2). 
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2. Nonperturbative Solutions of the Two-Body Dirac Equations 
For Two-Body Dirac equations of constraint dynamics applied to QED we have 



A{r)^--. (2.26) 
r 



For singlet positronium system we can obtain an exact solution [4] for the total CM energy w 



2 + 2/y 1 + ay{n + ^ (l + - - I - 1)2 

= 2m - ma^/in^ - ma'^/2n^{2l + 1) + ll/eima'^/n^ + Oia'^), (2.27) 

that agrees through order with standard spectrum found by perturbative treatment of the Darwin and spin- 
dependent terms in our Pauli form. Numerical triplet state calculations agree equally well with perturbative QED [5]. 

Many of the standard approaches to QED bound states have been applied in QCD in nonperturbative numerical 
calculations of the meson spectra without first testing them nonperturbatively in QED. Sommerer et al. [30] have 
shown that the Blankenbecler-Sugar equation and the Gross equations fail this test. This indicates danger in applying 
such three dimensional truncations of the Bethe Salpeter equation to quark models, for if failure occurs in their 
applications to QED how can similar non-perturbative (i.e. numerical) approaches based on the same truncations 
(but with QCD kernels) give results that arc trustworthy representations of the physics for meson spectroscopy? 



C. Two Body Dirac Equations for Meson Spectroscopy - The Adler-Piran Potential 

We obtain a constraint version of the naive quark model for mesons by employing a covariant adaptation of a static 
quark potential due to Adler and Piran [31]. From an effective non-linear field theory derived from QCD they obtain 

VAp{r) - A(C/(Ar) + Uq) {=A + S). (2.28) 

The original Vap is nonrelativistic, and appears in our equations in that limit as the sum of world vector and scalar 
potentials with 

AU{Ar < < 1) ^ 



In Ar 



VAp{r) = A[ciAr + C2 log(Ar) + ^ + ^ + 05], Ar > 2. (2.29) 

VAr Ar 



The explicit form for VAp{r) at all distances is given in [31] and [32]. 



1. Relativistic Naive Quark Model 

We reinterpret the static Vap covariantly by replacing the nonrelativistic r by ^/x^ = r, and parceling out the 
static potential Vap into the invariant functions A{r) and S{r) [7] as follows: 

A = exp{~Pr)[VAP --] + - + —, S = Vap + — ~ A. (2.30) 
r r r r 

(The constants ci, C2, C3, C4 are fixed by the Adler-Piran formalism while ei, 62 are the quark and anti-quark electric 
charges.) Thus at short distances the potential is strictly vector while at long distances the vector portion is strictly 
Coulombic with the confining portion at long distance (including subdominant portions) strictly scalar. Once A and 
S have been determined, so are all the accompanying spin-dependent interactions 

= $,((71,(72, Pi, A(r),^(r)); i ^ D,D' ,S0l,S02,SS,T,.. (2.31) 

Our bound state results are quite accurate, from the heaviest bottomonium states to the pion. They compare quite 
favorably with the results of Godfrey and Isgur [33] , but with only two parametric functions {A, S) as opposed to the 
six or so used in their approach. In the table below wc reproduce a portion of the entire spectrum given in [7]. The 
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quark masses and potential parameters are given by iriu ~ 55 MeV , rric ~ 1-5 GeV, ~ 58 MeV, A ~ 0.216 GeV, 
and AUo = 1.865 GeV. 



MESON 


EXP (GeV) 


(±MeV) 


THEORY 


— lie* 

r]^ : cc VSo 


2.980 


( 2.1) 


2.978 


tp : cc i oi 


o.uy ( 


(U.U ) 


Q 1 on 
o.lzy 


Xo • cc 1 ri 


o.OZD 


(r\ 9^ 

[V.Z ) 


o.ozu 


Xo • CC 1^ 1^0 








V 1 • 1 Pi 
Al • J X 


3.510 


fO 1) 


3.507 


Yn : cc l^Po 


3.556 


(0.1) 


3.549 


?7c : cc 2^50 


3.594 


(5.0) 


3.610 


V' : cc 2351 


3.686 


(0.1) 


3.688 


: cc l^Di 


3.770 


(2.5) 


3.808 


: cc 3351 


4.040 


(10.0) 


4.081 


: cc 23Di 


4.159 


(20.0) 


4.157 


: cc 33D1 


4.415 


(6.0) 


4.454 


TT : ud I^Sq 


0.140 


(0.0) 


0.144 


p : l^^i 


0.767 


(1.2) 


0.792 


bi : ud l^Pi 


1.231 


(10.0) 


1.392 


ao : ud I^Pq 


1.450 


(40.0) 


1.491 


ai : ud l^Pi 


1.230 


(40.0) 


1.568 


a2 •■ ud 1^P2 


1.318 


( 0-7) 


1.310 


TT : ud 2^So 


1.300 


(100.0) 


1.536 


p : ud 2^Si 


1.465 


(25.0) 


1.775 


7r2 : ud 1^Z?2 


1.670 


(20.0) 


1.870 


p : ud l^Di 


1.700 


(20.0) 


1.986 


p^ : ud l^Ds 


1.691 


(5.0) 


1.710 



2. Positronium and the Pion 

Positronium numerical spectral predictions of the constraint approach for hyperfine splittings are inadequate if we 
ignore coupling to the small (including the lower-lower one if)^) components of the wave function [5]. In the table 
below, Nc refers to the number of coupled equations, which for the fully coupled system is two for the singlet and four 
for the triplet states [5]. Units are in eV. As seen in the table, only the fully coupled system of equations (lower-lower 
and upper-upper for the singlet, the same in addition to tensor coupling for the triplet) produces accurate results to 
the require precision. 



I 


s 


3 


n 


Nc 


Perturbative (eV) 


Numerical (eV) 


Diff/^a" 











1 


1 


-6.8033256279 


-6.8032861579 


5.45E-02 











1 


2 


-6.8033256279 


-6.8033256719 


-6.08E-05 





1 


1 


1 


1 


-6.8028426132 


-6.8028074990 


-0.84E-02 





1 


1 


1 


2 


-6.8028426132 


-6.8028082195 


-4.75E-02 





1 


1 


1 


2 


-6.8028426132 


-6.8028239499 


-2.58E-02 





1 


1 


1 


4 


-6.8028426132 


-6.8028426636 


-6.97E-05 



The corresponding good n — p splitting obtained in [7] is spoiled if the we ignore these couplings, leading to 
^ 850 MeV, mp ~ 1060 MeV. The same relativistic structure in the constraint equations responsible for the 
success of the Sommerfcld-Balmer formula for positronium spin singlet states appears to be important for bringing 
the pion mass down to its observed value. 

3. Golstone Behavior of the Pion. 
As a bonus, we find [32], [6] ,[7] that the pion is a Goldstone boson in the sense that 

m^{mg ^ 0) ^ 0, (2.32) 
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while the p and excited tt have finite mass in this hmit. However, if the TBDE for the pion is truncated so that the 
couphng to the lower-lower component is dropped, then the pion loses its Goldstone boson behavior. Its mass no 
longer decreases toward zero with vanishing quark mass. This and the n — p result above support our contention 
that the pion docs not need to be treated in a special way insofar as the binding mechanism is concerned. The light 
pion mass as well as its Goldstone behavior is a natural outgrowth of the covariant Two-Body Dirac formalism. We 
now see how this model for the pion and other mesons holds up for a different probe, that of 2j decays. 

III. TWO GAMMA DECAY AMPLITUDES FOR POSITRONIUM AND QUARKONIUM 

Our treatment of decays in the sections below are for general angular momentum states but for illustrative purposes 
we begin by considering a treatment of singlet positronium or quarkonium systems. They can be viewed as bosons 
given by the state vector 



\iSo) = / d'pi;ip){bi'd%-bl'dl'^)\0). 



(3.1) 



Both the electron and positron (or quark and antiquark) are off shell but on energy shell. The amplitudes for the 
annihilation of a quark-antiquark pair into two photons are given by the Fcynman diagrams in Fig. 1. 



P+ s+- 

P- S- ■ 
P+ s+- 







P-- ki 





(a) 



P-- k2 



• ko e(«2) 



FIG. 1: Feynman diagrams for the annihilation of a quark-antiquak pair into two photons. 

The singlet amplitude for annihilation of a free e^e^ pair with momenta p-i- and p_ into two photons with polar- 
izations e"i,e"^ and momenta ki = (w/2, k), k2 = (w/2,— k) is 



M, 



(27r)3wV2 (p_ - fci)^ + - iO 

(P- — K2) + Tl^ — lO 

For positronium or quarkonium, we would replace this decay amplitude by 

' M 



(27r)'% 



(3.2) 



(3.3) 



Unlike free amplitudes, the fermion spinors and momenta in 15^,^27 S'l'e not on shell. 



A. Sixteen Component Two Gamma Decay Formalism 

The amplitude in Eq. (3.3) above is of the form (in CM) 



I d3pV(p)^[^(^+)(-p)r(p,k)u(-)(p)-«(-)(-p)r(p,k)u(^+)(p)] 

-^y'd3pVXp)Trr(p,k)[u(-)(p)«(^+)(-p) - w(^+)(p)z;(-)(-p)]. 



(3.4) 
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in which 



r(p, k) ^e^[7-e("-) ' T'jP ' 7-^'°-^ + 7-^^"-^ ' ^'^ + ''j T-^^"-^]- (3.5) 



(p — k)" + (p + k) + 

We replace this amphtude for general angular momentum states by 



d3pTrr(p,k)V(p), (3.6) 

where tpip) is our bound state wave function in matrix form in an arbitrary angular momentum state. Thus we 
are expanding our investigation from states to general ^Li and ^Li±i states. In the case of states what we 
are doing amounts to replacing the matrix wave function -0(p)[M^''~^(p)w^''+)(— p) — u'^'*+)(p)?;^'*-^(— p)] having a spin 
structure governed by free Dirac spinors by the matrix wave function ^p(Jp) which, unlike the solution constructed from 
the free spinors, is a solution of the full interacting set of Two-Body Dirac equations. Similar comments apply for 
the other angular momentum states. In terms of the Fourier transformed matrix wave function '4>{r) (defined below 
in Eqs. (3.27) and (3.31)) 

J d3prrr(p,k)V;(p) = J d'rTrmr) J d^p^^^Z^P-l) r(p, k)]. (3.7) 



Now 



/ 

exp(-^k.r) J d3,^^^^-^r,,(q)+ exp(*k • r) J d3p!^g-^rc(q), (3.8) 



,o exp(— ip • r) ^, , , 

' P (1)3/2 ^ r(p,k) 



where 



m-' + q 

Performing the Fourier transforms gives 



q = eM ^ „ , 3.9 

m-^ + q^ 

2 7-e("^)[m-q-7]7-e("i) 

^^(q) = e _2 2 ■ 



3 exp(-iq-r). 



(2^)3/2 rz.(q) 
=2. fl^.A^^)(r„, _ ,,^.^.exp(-mr) 



= e^^--ye'-°''>{m-i-fV)^^ '-f-e^"'', (3.10) 

and 

= e^Jly-,M(^ - ,-,-V) '"'P'-""-^ .e(°.). (3.11) 

This generalizes the configuration space form given in [15] to the amplitude below depending on the full 4x4 matrix 
wave function 

Mx^2, = J dVexp(-jk.r)rr[^>(r) J d'p^^^^L^2^r{p,k)] 

J dVrr{^(r)[exp(-ik-r)7-e("i)(m-i7-V)^^^^^t!!!!i^.e(«2) 

- exp(zk • r)7-e("^)(m - 17- V)^^^^^^^7-e("i)]}. (3.12) 
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The wave function will often display mild singularities at the origin typical for relativistic wave functions. For example 
the ground state solution corresponding to Eq. (2.27) is 

V'(r) = ^"'"^ " (r77za)+<i/^+V^/^'"') cxp(-amr/2). (3.13) 

y'47rr(2 + 2^/1/4 

The mild singularity at the origin appearing in this equation is rendered harmless by the smearing action of the 
Yukawa distribution that comes from folding the effects of the decay amplitude with that of the wave function [15]. 



1. Ax A Matrix Form of Solutions of the Two-Body Dirac Equations 



To accommodate the structure of the TBDE to the above decay amplitude we explicitly construct the 4x4 matrix 
wave function solution V'(r) of the equation. First we observe that one can write Eqs. (2.18) in terms of mass and 
energy potentials and their derivatives analogous to what is done in the case of two-spinless particles [5, 6], [34], [9]. 
In analogy to the solution (2.11) we gave to the third law condition in the spinless case we define 



Ml = TOi cosh L{S, A) +TO2sinhL(S', ^), 

M2 ^ m2 cosh L{S, A) + mi sinhL(S', A), (3.14) 

El = El coshg{A) - e2Smhg{A), 

E2 = £2 coshg{A)-ei sh\hg{A). (3.15) 

In terms of these functions the coupled Two-Body Dirac equations in an arbitrary frame have the form Siip = in 
which 

51 = exp [g] f3iEi ■ Vi + EiPij.i + M1751 - exp (g) '-^2 ■ d{gPi + i/?2)75i752, 

52 = - exp {g) P2^2 ■ V2 + ^^2/32752 + A/2752 + exp {g) ^Si • 9(^/32 + i^/3i)75i752, (3.16) 



with 



r, = p-^j:,-dgj:,. (3.17) 



The gamma matrices have block forms given in Appendix A. 

If we use the combinations (/)^ = V'l ± ^"4 ^-^d x± = '02 ^ V'sj then unlike Eqs. (2.24) and (2.25), the corresponding 
Schrodinger-like equations decouple [6] , [34] , [7] . We obtain [7] 

[p2 + 2m^S + 3"^ + 2e^,A - A^ 

^^v^g + ^g'^ -{g' + L'f + g'F' 

_ L ■ {ai + <72) ^, ^ ^ . _ ^ ^^/^ . ^ 
r 

-\-2F'if ■ p + iK'{(Ji ■ f(j2 ■ p + (72 ■ r<Ti ■ p) 

+ai . a2{\^^g + ^L' - \g"^ - g'F') + ai ■ fa2 ■ f{^W^L - (|- + F'))]4>+ 
= b^w)^+, (3.18) 
where the prime symbol stands for d/ dr. We have used the abbreviations 

F ^ ilogp-g, 
V = E2M1+E1M2, 

{g + L) 



K = 



2 



l\r) = _li^.M2-i^.M, 
^ ^ 2r E2M1 + E1M2 

1 E1M2 - £^2^1 , , . 

2r E2M1 + E1M2 
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We work in the CM frame in which P = (1,0) and r = (0, f ). Once we find the four component solutions (p^ to this 
equation we can obtain the other twelve components x_(_. In Appendix B we find from Eq. (3.16) 

X+ = ^^{Af2[<Tl-p-^<T2-V(-g-L + g^l.<T2)] 

~Ah[cT2 • p - • Vi-g -L + ga^- cr2)]}(^+, (3.20) 



and similarly 



and 



in which 



and 



i 



-Ei[(T2 • p - -<Ti • V{-g -L + gtTi- cr2)]}0+, (3.21) 



iE2Ei+M2Mi) ^ 

V 

-^[{E2D^+ - E,D2+)^{M2Dt+ - MiDt+) 



(A'hD^- + M,D^-)l{E2D++ + E,D++)]^+, (3.22) 



D++ = exp{g)[arp-^cT2-V{-g - L + g<Ti-a2)], 

= CXp{g)[(Ti-p+-(T2-V{g - L - g(Ti-(T2)], 

D^- = cxp{g)[(Tvp+^(T2-V{g + L - g(Ti-(72)], (3.23a) 



D++ = cMG)[(r2-p-^(TrVig - L + g(Tra2)], 
D^^ = exp{g)[cT2 ■ p + -a-i ■ Vig - L - gai ■ (T2)], 

D^- = CXp{g)[(72-p+^(Ti-V{g+L-g(Ti-(72)]. (3.23b) 

We then further define four component wave functions ip^. , 77^ related to the above by [9] 
(f)^ = exp{F + Kcri-r(T2-r)tp^ = exp F{coshK + smh.K(Ti-r(T2'^)'4'±j 

X± = exp{F + Kcri-r(T2-r)ri^ = expF{coshK + smhK(Ti-r(T2-r)ri^, (3.24) 

In this case the decoupled form of the Schrodinger-like equation for ?/'_|_ has the convenient property that the coefficients 
of the first order relative momentum terms f • p and (cri-f(T2-p+<T2-fcri-p) as appear in Eq. (3.18) vanish. We obtain 
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[9] 



in which 



, , , , 2F'(cosh2i4'- 1) „ „ 2i5:' sinh 2ii: 

r r 



L • {(Ti+(T2)[- 



F' F'(cosh2X-l) (cosh2A"-l) K' smh2K 



r r j.^ ^ 

+ L • (cTi -0-2) (/' cosh 2K - q' sinh 2K) 

+ iL-(TiX(T2 ((?' cosh 2K + I' sinh 27^) 

„ , ii:'(cosh2A'- 1) sinh2ii' K' F' smh2K , 
+ cri-rcr2-rL • (o-i+(T2)( ^ H — H ) 

„ , , F' sinh 2 A" F' (cosh 2 A' - 1 ) A'' sinh 2 A' 

+ <Ty(T2[k{r) ^ + 

Zr r r 

_^ ii:'(cosh2A:- 1) ^ sinh2A' (cosh2is: - 1)^ 

_ , , 3F' sinh 2 A' F'(cosh2X - 1) , , A'' sinh 2 AT 

-(Ti-f(T2-f [n(r) + \ i + 2F'K' 

r r r 

■iK' {cos\i2K - I) . 3sinh2A: (cosh 2A" - 1) , 

^ — - ^ + — -2 — + ^ — -2 — 

= (3.25) 



k{r) - 2V g--g --GK -27 + -, 

n(r) = \/^K -^V^g -2K'F' + g'F' - 

m{r) = -l\/^g+^g'^ + g'F' -K'^, (3.26) 

For equal mass singlet states, the hyperbolic terms cancel. The spin-orbit difference terms in general produce spin 
mixing. 



2. Matrix Form of the Wave Functions 



We now construct the 4x4 matrix forms of the wave functions (appropriate for a spin-one-half particle-antiparticle 
system) from the sixteen component forms (appropriate for system of two spin-one- half particles). We begin by 
writing the 16 component spinor wave function as 




(3.27) 



The spinors ijjj^ as well as (/)_|_ = Vi ^ "04, X± = "02 ^ V'3 £^re themselves four component Pauli spinors (upon 
which aii,a2i operate). The conversion from sixteen component spinor wave functions to four by four matrix wave 
functions now can be carried out in a two-step process. First, as in [34, 35], the "energy" or q space column vector 
direct products are converted to 4x4 matrices as follows (recall the factor of iay plus the transpose operation changes 
particle spinor into antiparticle spinor) 



(1) 



(2) 



*(l)\E'f2)«"2 = *(i)*f2)'?l ® «Cr2, 



(3.28) 
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in which aQ,ai,qo,qi,;i = 1,2,3 are the 2x2 unit and three Pauh matrices in commuting spaces (spin and energy 
space) 

= StjQo +i£i]kqk, (3.29) 

and whose direct products form the Dirac matrices. Second, the (l)±,X± four component Pauh spinors are converted 
to 2 X 2 matrices in a and q space by 

X± X± = {x±o<^o + X±-o-) ■ (3.30) 
Together, the 4x4 matrix wave function in cr, q space is 

(^ (8) lo+'Y' ®93 + ^®gi + ^® iq2)qi ® i<J2 

= (^ (8X71 + ^ ® + ^ ® go + ^ <8) g3)«cr2 1 

( ^ '^gi+ ^ '^»g2+ ^ '8'go+ ^ t^gs) 

^ (^®'7i + ^®«g2 + ^«)go + ^®93), (3.31) 

where for convenience we have absorbed the factor icr2 into the wave functions as the wave function is arbitrary up 
to a constant multiphcative matrix and we have used the same symbol for each of the transformed wave functions to 
simphfy notation. In our work below we drop the direct product symbol it being understood to apply whenever 
a and q space matrices multiply one another. 

The four component spinors ijj^. and rj^ are similarly transformed into matrices which can be expanded in terms 
of (Tq s-nd cr. 

'n± V± = ('?±o^o + i7±-o-) ■ (3.32) 
With A a generic matrix, Eq. (3.28) leads to 

o-i • AtA^ A-cr (tA+oCTo + V'+'Cr) 

= A • '0+o'o + A-cr?/)^o + * A X 

(T2-A'4>_^_ -A ■ [ip_^_Qao + xp^-cr) a 

= -A ■ xp^ao - A-crijj+o + iA X ip^-cr, 

(Ti-A(T2-Aip^ -cr-A (V'+QfTo + V'+'Cr) cr-A 

= -A^^+oO-o + (AV+-2A-t/>+A)-(T, (3.33) 

which are needed to convert Eqs. (3.20), (3.21), and (3.22) into their matrix counterparts. In terms of matrix wave 
functions 0^, x±, V'±j ry^, we find that Eq. (3.24 ) becomes 

0±o = exp(F-if)V^±o; 0± exp(i^ + X)(l - (1 - exp(-2i^))ff ) • 

X±o = exp(F-if)r;±o; x± - cxp(F + if )(1 - (1 - exp(-27^))ff ).»7±. (3.34) 

We write the 4x4 matrix wave function form -ijj of the sixteen component tp to be used in our decay amplitude in 
terms of the matrix form 

*(r)=^^(V'+'7i + i'-iq2 + ??+go+'7-93), (3.35) 

where ip^ — ip_^_QaQ + tp^-cr is the 2x2 matrix form of the solution of the above Schrodinger-like Pauli equation 
(3.25). 
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Using these four components, the remaimng twelve components ip^a, V'-i ''1±oj ^^'^ 'n± ^^'c obtained from Eqs. (3.20), 
(3.21), (3.22), and (3.34). In aU of our decays the particle antiparticle pairs have the same mass: nii = mi = m 
and so £i = £2 = e = w/2. Using the definition 

Ml = A/2 = A/ = mexp(L), 

El = E2 = E = seyip{-g), (3.36) 



we show in Appendix B 



V+o = ^^P(g^+ 2-^) [p _ ^ v(L + 2F + 2/0] ■ [1 + Q,„ff] ■ 

r7+ = ^^^^^{(P - ^ VL) + Q,(ff . p - ^VL)}V'+o, (3-37) 

and 

V-o = 0, 

V = -^^^[1 + Qpff] • [iP + 1^{L - 20)] X [1 + Q„ff] • V+. (3.38) 
The final four components of the four by four matrix wave function found in Appendix B are 



° ^ 2EM 2ME ^ 2 ^ 2 J-^^+O' 



+ A/2 , exp 25 , 
^" 2i;M ^+ 2EM L ^"^P J 

[p - ^V(L + 65)] [p - ^V(3L - 25)] X [1 + Q,„ff] • 

p - ^V(i + 25)] X {[p - '-V{L - 25)] X [1 + Q™ff] • V+} ) , (3.39) 



where 



Qp = exp(2/C)-l, 

= exp(-2/C)-l. (3.40) 

We also show in Appendix B how for both singlet and triplet states these solutions together with Eq. (3.35) are 
related to the solutions governed by the free Dirac spinors in the absence of interactions (see also Eq. (3.4) and 
discussion below (3.6)). 



3. Covariant Normalization Conditions for the Matrix Wave Function 



In this section we discuss how the norm of our matrix wave function will differ from the naive form of 

i j d^xTrg^-^^^^^ J d^xTr^i^l^+ + ^pltp_ + rjlij^ + rjlrj_) = 1. (3.41) 

In a series of papers in the context of constraint dynamics, H. Sazdjian has shown [36] how this norm must be modified 
so that, like its nonrelativistic counterpart, its constancy is connected to a conserved, in this two-body case, tensor 
current. The norm he developed was not for the solution of a quasipotential equation like Eq. (3.25) but rather 
developed from a set of two-body Dirac equations similar to those we use here. It deviated from one like the above 
by terms that depend on the interaction as well as the way in which the interaction depends on the CM energy. Later 
work [37] simplified the norm to one that is interaction independent when the interaction is independent of the energy. 
In terms of the 16 component spinor solutions of the Two-Body Dirac equations given in Eqs. (2.18) we found the 
norm condition of 

d^xlijUl + 4u;2/3 /J )^] = / d^xtfj^'Ctp = 1. (3.42) 
ow^ J 
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If the matrix A is CM energy independent, then the norm is like that of the (one-body) Dirac equation (with no 
energy dependence of the interactions). We caU the norm of Eq. (3.41) the naive norm (NN) and that of Eq. (3.42) 
the two-body Dirac norm (TBDN). The connection between the matrix interaction function A and the core scalar 
and vector interactions appearing in Eqs. (2.18) were found in [37]. There we showed that 2.18 has the hyperbolic 
structure 

Sii: = (cosh(A)Si + sinh(A)S2)V' = 0, (3.43) 
S2i' = (cosh(A)S2 + sinh(A)Si)?A = 0, 

in which 

SiV' = (5io cosh(A) + 520 sinh(A))^/' 0, 

82^- = (^20 cosh(A) + 5io sinh(A))^/' ^ 0, (3.44) 

with 

5io-0 ( - Pi^i ■ P + ei/5i75i + '^i75i)V' 

^20^ = (/32S2 • p + £2/32752 + ^2752)^ (3-45) 

and 

A = ^751752 [i(a^i) - 7i ■ l2Gixi-)]- (3.46) 

with L and Q given in Eq. (3.36) (see also [41]). In matrix form the connection given in Eqs. (3.34), (3.31), and 
(3.27) between the matrix form of the wave function t/j of (2.18,3.43) and 'I' is 

■0 = exp(F)[cosh A'*(r) - sinhA'S • f*(r)S • f 

= /C*(r) (3.47) 

In Appendix C we show that in terms of the matrix wave function 4' solution (3.35) to Eq. (3.25) the nomalization 
condition can be written as 



J SxTrijj'^Cip = j d3a;rr(/C*(r))'^/:/C*(r) = 1. (3.. 



There we also give the matrix form of the operator C. The deviation of the matrices K, and C from the unit matrix 
will affect the decay rates. The decay amplitude (3.12) in terms of the matrix wave function 'l'(r) is 



Mx^2^ = e^yi j (frTr„,{lC^>{v) 

X [exp(— ik • r)(j3(7icr-e'"^-' (to — iq^qicr-'V) ^ ^-q^qicr-e^"^^ 

r 

+ exp{ik ■ r)g3(/icr.e("^)(TO - zgagicr- V) '^"''^^^""'"^ gsgiT-e^"^)]}. (3.49) 



4-. Scalar and Vector Wave Functions in Vector Spherical Harmonics 

Given the above wave functions we now write down the total 4x4 matrix wave function in terms of -0+. The spin-zero 
part of the total wave function is governed by '4'+o^ spin-one portion by xp_^_. These wave functions appear in the 
forms 

w, - iMv- 

rP+ = %i^Y,™+ + %i^Y,„,_ + ^X,„„ (3.50) 
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where the labels on the radial wave function u refer to the Isj quantum numbers of the solutions to Eq. (3.25) and 

Yjm+ = {a+r+rb+p)Y,jm, 
Yjm- = (a_f+r6_p)Yj„,, 

^jni = / ., ^ , -,V ' (3.51) 
VJU + 1) 

arc vector spherical harmonic cigenfunctions of with eigenvalue l{l + \) where I = j + 1, j ~ \, j respectively. The 
coefficients are 



J + 1 / J 



2j + l' V2J + 1' 



In our work below, there will be no spin mixing and the unnatural parity solutions {u^-y-/r)lLjm{^) will not contribute. 

For spin-singlet states [xpj^ = 0), Eqs. (3.35), (3.37), (3.39) imply the following combination of scalar and vector 
wave functions 

^'|s=o = + ip-Qcrniq2 + V+ ■ <^<lo)- (3.53) 

In Appendix C we show that the TBDN for spin-singlet states is 

1 = \ y"d3xexp(2F)([exp(-2A')(^^oV'+o + ^LoV^-o)+cxp(2/Or7t_-r,+ -2sinh2A'r,^-frj+-f] 
dL 

-f [2i«2_[cxp(-2i^)(V'%V'+o - ^lo^-o) - cxp(2A0r,; • r,+ + 2sinh2A-r,^ • vr)+ ■ f] 

+27i;2|^[2exp(-2i^)(2^>;o^+o + V^t_o7^_o)]) (3.54) 

The naive norm (NN) is given by the above two-body Dirac norm by £, /C ^ 1 or equivlently by exp(_F), exp(A') — >■ 1, 
dL/duP', and dG /dw^^Q and is 

i J dMl^UV'+o + vl-V++ i^U^-o] = 1- (3.55) 

Appendix D gives from Eqs. (3.37) and(3.39) the needed radial forms for the contributing wave functions in terms 
of the radial portions of the solution to Eq. (3.25). It requires the radial form of Eq. (3.25) which is simply 



r dr^ r2 2 ' ' r 



;r + + ±v^A - ^ [VLf]^ = exp(-2g)]^, (3.56) 

where [411 



22 „,,„/- o/7^ _ C I c2 I /I /l2 



-6''exp(-2g) = 2m^S + + 2e^,A- A\ (3.57) 
Appendix D gives us the relations between the contributing wave functions V'+oi ^-O' They are 



and 



/+ — \ ^ ^ J7n— I ^ -"-jm+y? 



£■ ' ' dr 2' r ' r V 2j + 1 
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For the norm wc also need 



For spin-triplet states ("0+0 ^ 0)' ^Iqs. (3.35), (3.37), (3.38), and (3.39) imply the combination 

= ^^(■0+-o-gi + ip_-(Tiq2 + jy+o'Jo + V-'f^Qs), (3-61) 

and the contributing wave functions are xp_^^,xp_, V+o^ ^^"^ V ■ Appendix C we show that the TBDN for spin-triplet 
states is 

i J d^xcxp{2F){[cxp{2K){^(,l■^P+ + + ril q_) + cxp(~2A')r/%??+o 

— 2 sinh 2/v (■0^-f'0_|_-f -|- ■(/'- ■f''/'- ■ r+rj^^ rri^ • f )] 
dL 

+{[2u;2_[exp(2A')(-0t-^+ - ^pl■^|J^ + rjl-V-) ~ exp(-2X)7?%77+o 
— 2 sinh 2A'(i/'^-f'0+-f — i/j^-fi/)- • r+r]^ rr]^ ■ f)] 
+4^«'^([-exp(2/0(V'^V'- + »7--»7-) + 2cxp(-2if)ry^or7+o 

+2smh2K[{xp\_-rip+-r + ipl-rip^-r + r]^_^-rr]^ ■ r)]]}) =1, (3.62) 
while the naive norm (NN) is 

i J d^xixPl-'iP+ + + 7;%j;+o + »7--»7-) = 1- (3-63) 

In Appendix D we show that Eq. (3.39) gives ■0_ from 

in which 

^ = T^f{<&---2(j + l)62exp(-2g)+2V7(7TT)a>+_ + ^ + ^ + C„„„ + :^^ 

2j + 1 \ J, r dr dr r 

r dr dr r 

- 2jB' exp(-2g) + 2V7(7TT)a>_+ + % + + C,,+^^ + Gpp-f }^Y„„+ 

r r r dr dr r 

+{-*+- + V7(JTT)[2<i>__ - 262cxp(-2C;) + ^ + ^ + C,,„ + + G,„,-^]}— Y,„+) , (3.65) 

r r dr dr r J 

(see Appendix D for explicit forms of the functions Amm, .., Gpm,). This equation requires the coupled radial wave 
equations for spin triplet states that follow from Eq. (3.25). They have the form [9] 

r dr j,^ J. J. J, 

[_i^, + Z(l_l)+ci,__]!l(2^ + <j,_^!l(W = s2exp(-2C?)^^S^. (3.66) 
r dr'^ j,z J, J. J, 

(See Appendix D for the explicit forms of <i>±±). Thus with Eq. (3.50) we have 



+ V7(7TT)[2a>++ - 2S2exp(-2g) + llfi + + a„p + + G„,pf ]}^Y,„,_ 



^ M«J_i)ij exp(2g) 



E r 2EAI{2j + l) 



- 2(j + l)62exp(~2g)+2V7(jTT)*4 



5 1 r Omm + ( r LfmmJ^I 

r r dr r 



-[<!>-+ + v/KJ+1)(2$++ - cxp(-2C;)+^ + ^ + C™P + + G™p)^)]^^}, (3.67) 
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and 

^U+m3_^^^U+m exp(25) 



E r 2EM{2j + l) 

dr r 



{[-$++ - 2jB^ exp(-25) + 2v/j(.? + l)*-H 



+ y7(JTT)(2<i>__ - 262cxp(-2g)+% + ^ + + + Gp™)-^)]^^}. (3.68) 



In Appendix D wc also show that 



and 



r dr r 



V+o = ^^^Yim. (3.69) 



q = x^Xj-™, (3.70) 



in which 



^ = ^xp(g + 2K) ij - 1) - 2Q^ _ ^^^^ ^ d ^ 5 ^ r7""^-^i)i.- 

r h r dr 2 \i 2j + 1 r 



and 



V 



+ 



CXp(g) d j(Qm + 1) - 1 , J + 1 



r 



+ + r (^^+ 2 ) + l r ' ^^'^^^ 

We use these wave functions to compute composite 27 decay amphtude Eq.(3.49) which after performing the q 
space trace gives 



Mx^2^ = j d^rcxp{F)Tr„{cxp{^ik ■ r){z[cosh/i - sinh Xcr ■ fV'.cr ■ f]cr ■ e'-"^\cr-V f^'^^^ "^^-^ 

+ [cosh A'77_^ - sinh A'cr • ny+cr • f]cr • e^^^^m ^^'P^"'"'''' }cr • e^"^) 

+ exp(ik • r){i[coshiir-i/)_ — sinh A'cr • Y%jj_cr ■ f]cr • e^"^-'(<T-V) 

r 

+ [coshA"ry_^ - sinhAV • vr^^a ■ r](T ■ e("=)m^^^^^^-^}cr • e'^i)). (3.73) 
The trace ehminates the contribution of the portion ijj^{v)qi of the wave function. 

5. Decay Amplitude for ^Li Composites 

Substituting Eq.(3.53) into (3.73), 

MiL,^2-f = -^-^ J dVexp(F)Tr4{zexp(~AOV-o 

x[exp(-ik • r)cr-e("i)(cr-V)cr-e("^' + exp(ik • r)(T-e^"'-\(T-V)(T-e^°''^]} 
+{[exp(A')j7+-cr — sinh(Ar)2f -rj+cr-f] 

x[exp(-ik • r)cr-e(°i)mcr-e + exp(tk • r)cr-e^"^)mcr-e("i)]}) '^'"P^"™'^) ^ (3.74) 
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performing the remaining trace gives using (3.58,3.59) gives 

MiLi~.2-f = -^e^ e("i) X e^"^) ■ ^ dVexp(F)[exp(-ik • r)- exp(+ik • r)] 



.^^^^^,K,r,..K-''i^ (3.75) 

with unit vectors defined in terms of the photon decay momenta and transverse polarization vectors 

z= k, 

e^±). (3.76) 



(x±iy) _ (±) 



The integral forms appearing in Eq. (3.75) are treated in Appendix E in which we show (with g(r) appropriately 
defined) 

y dVexp(-ik-r)f5(r)Y,„(fJ) 



4^ (-^)/^(2L±l)(ji;00|/0)y° 



x[k(jl; 00|/0)(5,„o - e^^^jl; -ll|/0),5™_i + e(+)(jl; 1 - ll|/0)d;„i]. (3.77) 



We also show (with /±(r) appropriately defined) 
y dVexp(-2k-r)/±(r)Y,>„±(r!) 



J2 i-^y' O'l; 001/0) / drr^f±{r)[{a± ~ 2ib±)j,, (kr) - ih±j'y {kr)kr] 
x[k(jl;00|/0)<5„o - e-(jl;-ll|j'0)<5„_i + e+(jl; 1 - ll|/0)<5™i], (3.78) 



in which 



j^ikr)^ ^^^-^^''^-^^_^^^^-^^''\ (3.79) 

Thus 

MiL,^2-, = ~y27TTe("^' X e("^) • k{Fj=,(l + (-)^),5™o (3.80) 
+ E ^?=1(1 - (-r V^ ' 00|/0)(jl; 00|/0)J™o}, 
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in which 

Fj^i = — 2«7re^(— i)-' / drmrexp{~-mr)jj{kr)krcxp{F + K) 



S + iy 2.7 + 1 r J^2j + 1 r ' 

poo — 

of J = -27re2(-i)j' / drexp(-mr)exp(F)(jjv(fcr){(mr + l)exp(-/i:)^ 



-2mr shMK)i ^ + , ^^)} 



+ [(1 + ^)jAkr) + ^fAkr)kr]^ ^-l^'^^^^f^jeMK)). (3.81) 

Notice that this amphtudc (3.80) is zero for j odd consistent with the Landau- Yang Theorem. We caU this amphtude 
the two-body Dirac amphtude (TBDA). What we caU the naive amphtudes (NA) would correspond the use of the 
naive norm (/C = £ = 1) together with exp(F), exp(iir) ^ 1 in Eq. (3.81). 



6. Decay Amplitude for '^Li±i Composites 
Using Eq. (3.61) in (3.73) leaves us with 

A^3i^^,^j_,27 = -y V^i^ / (i^rexp(F)Tra^(cxp(-ik • r){[iexp A'i/)_-cr 

-2sinhA'Vj_-fcr • f]cr • e^^^'cr-V + mcxp{-K)r]^o(T ■ e("i V • e'"^)} 
+ cxp(+ik • r){i[exp Kxp_-cr — 2sinh A't/j-'fcr • fjcr • e^^^-'cr-V 

+mexp(-A')7] ,0*^ • e^"'^}*^ ' e^^^)}) ^^^^i^^. (3.82) 

r 

Notice that only two of the four portions of the triplet wave function (3.61) survive that trace. Performing the a 
space trace and using Eqs. (3.64) and (3.69) together with 



we obtain 

'exp(— ?7ir) \ ' 



J (i'^r[cxp(— ik • r)+ exp(zk • r)] cxp(F) ^ 



_^ "(j~i)iJ [Y^.^^^_(^) .g(aOf.g(a2) ^Yj,„_(12) -e^^^^f-e^"^) - Yj>„_(») -fe^^^^-e^"^)]} 



. + [„exp(-A-)^ - 2(m + 1/,-) sinh A-( Jiii^it™ + 5c«il7_2_)y,„.,t2)e(«.) - ,i«>BEt!^' 



With Eq. (3.76), the shiiplest terms in the above expression include forms hke 

dVexp(-zkT)5(r)y„„(n) = 47r(-^)■'y„„(Ofc) / drr^jjikr)g{r) 



^ V4^(2j + l)(-z)^5™o / drr''j,ikr)g{r). 

Jo 

Stepping up in complexity we have the transverse parts of the dyad form 

J d^r exp(— ik • r)G(r)f f . 
In Appendix E we show that transverse portion is 

(l-kk)-^ d^rexp{-ik-r)G{r)rr- (1 - kk) 

-J Jo 



J' = li-2| 

+(,eH..H,,y?^eZ±i) g " + ';"'"'' (,2;,0|/0)02;-22|/,)..„_. 

i' = li-2| 

(j2;00|/0)0-2;2-2|/0)(5™2} 



j' = \3-2\ 

/•oo 

X / r^dr{-iy' jy{kr)G{r). 



Finally, wc need the trace as well as transverse parts of the dyad forms 

J dVexp(-^k•r)F±(r)Y,>„±(^7)f 
= j d^rexp{-ik-r)[a±F±{r)rY^,„+b±rF±{r)pYj^)r. 
In Appendix E we show that with Eq. (3.76) the trace portion of Eq. (3.88) is 

(eWe(-) +e(-)e(+)) J exp(-ik • r)F±(r)Yj„,± (f^) • f 



)a±^ATr{2j + l)i-iySrnO / r^drF±{r)jj{kr) 
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while the transverse part is 

(1 - kk)- j dVexp(-ik • r)F±{r)Yj,n±{^)r ■ (1 - kk) 
1 f°° 

= (e'+'e'-) + e(^)eW)-v'4^(2j + l){J„o / r^dr{-iy j,{kr)[{a± + 3tb±)F±{r) + i6±rF4(r)] 
•J JO 

V ^ii^-i^^O-2;00|j'0)0-2;00|.7'0)<5„,o r2dr(-z/j,v(fcr)[a±F±(r) + i5±rF4(r)]} 



J+2 

■ E 



{„H,,,.>,;851l±i) f " + 001/0)02, -221/0)...... 



j'=li-2| 

J+2 



+ „.->..-),/E?Z±i) ^ (l±(_l£i)(,,;00|/0)(,2;2-2|/0)..„,} 

nOO 

X / r2dr(-i/jjv(A:r)[a±F±(r) + i5±rF4(r)]. (3.90) 
Jo 

After integrations by parts, substitution of vahies of a±,ib± and combining with the other portions, we obtain 



+ y(2jTT) 0-2; 001/0) + {-y')2Bff^Je^"^^-e^'^^\j2;00\j'0)d^o 

j'=b-2| ^ 

-V6[e("^) • (e(+)e(+)) • e*"^' (j2; -22|j'0)<5™_2 



in which 



. (e(-)e(-)) . e("^)(j2; 2 - 2|/0)J™2]) , 



Aj=i±i = — - — f / drcxp{~mr)cxp{F)(^{—3jj{kr){mrexp{—K)—^ 



(3.91) 



2 



+ (mr + 1) cxp_R'{[(/ (fc7') H j'Jkr)kr) 

7 + 1 



j + 1 V 2i + 1 r 



+ (-,i(frr) + j/(t'-)tr)l/ ,/^j°";"'-' l) 



= f / dr exp(— mr) exp(F)(mr + 1) 

3 Jo 



:(expif{[(-A__i)j.,(fc^) + _!_/., (fcr)fcri + ^ 



2 j + 1 r 



+ [(- + l)j,,(fcr) + if ikr)kr].U—^^^^^} (3.92) 
J J V 2.7 + 1 ?■ 



As in the singlet case we obtain zero amplitude (3.91) for odd j. We also call these amplitudes the two-body Dirac 
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amplitudes . Again, the corresponding naive amplitudes would correspond the use of the naive norm (/C = £ = 1) 
together with exp(F), exp(X) ^ 1 in Eq. (3.92). 

B. Decay Rates 

From the above two sets of amplitudes we construct the decay rates. In our present case, we have 

w 

w 

b = Ip^I = -. (3.93) 

Also, we are not interested in the decay of a state with a definite magnetic quantum number. Rather we are interested 
in the average over all m. The Lagrangian that leads to the Feynman amplitude for the decay process is Lorentz 
invariant. Consequently the amplitude and our bound state adaptation conserves total j, m. This implies that we 
can sum over final states in an unrestricted way that is most convenient, without picking only special helicities that 
one expects to contribute. The details of the amplitude should do this automatically. Using the general decay rate 
formula [38] we obtain 



e(°l).e(<»2) 



1 



{2j + l)16(27r) 



dQk \^x^2^\\ (3.94) 



E(c»l),e(a2) 



in which we carry out the initial state m average and final state polarization sum independently. For spin singlet 
states with amplitude (3.80) this becomes 

^(^^'-^^^- (2, + l)W)J -^"- E k'-^x.("^).k|^ 



J + 1 



x(2j + im^iil + (-)^>™o + - (-)-'')^^^^^^(il;00|/0)(jI;00|/0)5„oP 

^^\F,=iil + i-y)Srr.o+ E ggkl - i-y') ^^ ~ \ jl;00\fO){jl;00\fO)S„,o\'- (3.95) 



We have summed over the following four independent polarization combinations 

. k = e("2) . £. ^0, (3.96) 



with 

Y |e("^^ X e("^) -kl^ = J2 [1- |e^"'^ • (e^"'^)*n =2. (3.97) 

e(<»l),eC«2) e(°l),£(°2) 

Note that only the zero helicity states (corresponding to both photons being either left or right handed polarized) 

e("i),e("^) = -L(x + ,5^),_L(x_,y)^e(+),e(-), 

e("^),e("=) = -L(x-^y),i=(x + ^y)^e(-),6(+), (3.98) 
give non-zero contributions to the rate factor 1 - |e("i) • (e("2))*|2^ xhe total helicity ±2 states 

e("^),e("^) = -^(x + iy),-^(x + iy)EEe(+),e(+), 

e("^),e("^) = -^(x-zy),-^(x + zy)^e(+),e(+), (3.99) 
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give zero contribution. Performing the angular integration gives [42] 



rCSo ^ 27) = j^\Fo + G^^^(01;00|10)|2 = j^\Fo + G^^^\^ 



(3.100) 



and 



rCD2 



2j) = -^\2F, + ^ 2G;ii(l-(-)-'-')0-l;00|/0)0-l;00|/0)|^ 
^ |i^2 + G^^^(21;00|10)2 + G^^^(21;00|30)2|2 



(27r)4 
1 



1^2 + -G^- 



(1) 



(3) 1 2 



Using Eq. (3.91) for triplet states ^Li±i our rate formula is 



(3.101) 



r(^Lj=!±i 



27) = ^ E 

^ ^ m,e(°i),e(°2) 

J+2 



A,=,±ie("i) •e("-)(l + (-)^)5™o 



E + 001/0) 



j'=b-2| 



X (^i?]i;±i2[6("i).e("^)(j2; 00\fO)S^o 
-V6[e("i) . . e("=)(j2; -22|j'0)5™_ 



+e("i).(e(-)e(-)).e("2)^j-2.2_2|/0)<5„ 



(3.102) 



Notice from Eqs. (3.98) and (3.99) that this rate in general includes both helicity zero and helicity two contributions. 
In the case of '^Pq decay we have j = m = and so performing the polarization sum gives [42] 



r^Po ^ 27) = 
1 



J2 |Ao + 2i?(^)(02;00l20)2l2l6(-).e 



(Q2)|2 



(2,),|Ao + 25, 

This rate includes only helicity zero contributions. 
In the case of '^P2 decay we have [42] 



(2) ,2 



(3.103) 



r(3p2 ^ 27) = 



1 



2(27r)4 



E 



mO 



m,£(°i),e<°2) 

+ (22;00lj'0)2B^^"'('e("i'-e("^)(22;00lj'0)(5 

-V6[e^"'^ ■ (e(+'eW) • (22; -22l/0)(5„ 
. (e(-)e(-)) . e("^)(22; 2 - 2lj'0)<5™2]] 



(27r)4 



\A2 - 2[B!f^ (22; OOlOO)^ + S^^^ (22; 00120)2 + B^^^ (22; 00140)2] p 

+12] [B^°^ (22; 00)00) (22; -22)00) + ^ (22; 00)20) (22; -22)20) + b!^^ (22; 00)40) (22; -22)40)] )2 
+ 12) [S^"^ (22; 00)00) (22; 2 - 2)00) + ^ (22; 00)20) (22; 2 - 2)20) + B^"^^ (22; 00)40) (22; 2 - 2)40)] )^ 



i2ny 



9r(0) r,p,(2) ocr(4) o(0) 
|^2 + ^^ + ^^ + _^|2+24) 2 2 



(2) 



3B. 



(4) 



7 



35 



7 



35 



(3.104) 
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1. Positronium Decays 

For these decays we ignore the effects of the potentials on the norms and ampUtudcs since they are relatively weak 
(/C = £ = 1 or exp(i^), exp(/-i:) ^ 1). 

a. ^Sq Decay The amplitude for ^5*0 positronium decay is from Eq. (3.80) 

Tis,^{Fo + gI'^), (3.105) 

where for the weak potentials we expect in QED (with exp(F), exp(if ) 1) 

f°° u+ (r) 

Fq = — 2i7re / drmr cxp{~mr)jo{kr)kr—^ (3.106) 

Jo ^ 

G^^^ = i2TTe^ [ drexp(-mr)((mr + l)[ji(fcr) "°°"^''^ + mr[ji{kr) - j'o{kr)kr f^^°^''^ 
Jo \ "f" ^ 



where we use 



r E r E r w^l + 2a/ [wr) r 

and 



7+0 _ cxp(g) d L' u+„^ _ exp(g) d M+o _ 2 d u^^q 



r E dr 2 r E dr r w {1 + 2a / (wr)) dr r 

This wave function is one of the small component ones. For positronium, w — 2ra + 0{a'^) and so 

— = J 1 + O a )) 

and (with k = m(l + ©(a^)) 



Fo = 


POD 

—2i'Ke? / drr^m 
Jo 


(0) _ 
— 




— 27re^ / o?rexp( 
Jo 



exp(-mr) r d 

Ji[mr)mr{ ■ — )— lAooo 

mr + a dr 



f drexp(-TOr)((mr + l)[ji(?7ir), / Vooo 



m dr 

with the nonrelativistic wave function given by 



-'000/ 

2 / 1 V 2 



1 

2 

In Appendix F we obtain the well known form for the decay rate: 

G, 



'"o 



(27r)4 m2 2 

and show that the small component portion Fq does not contribute to the singlet decay rate at this order. 



(3.107) 



,^ -^:2m_ (3.108) 



+mr[ji{kr) + j'^{kr)kr]-^%^^) (3.109) 



(ma)3/2 i?(r) 
V'ooo = exp(-amr) = — = (3.110) 



replacing the relativistic one u^QQ/r. The NN (3.41) becomes 

drr^[('^\ +(^] (3.111) 



-1^(0)1^^ = ^. (3.112) 
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b. '^Po,2 Decay The branching ratio for these decay have not been measured since the decays of those states is 
so largely dominated by the dipolc transition to the '^S'l state. Nevertheless, it will be of value to determine if our 
covariant formalism yields the standard results given in [39] and [40]. The rclcvent amplitudes given in Eq. (3.92) for 
weak potentials (£,/C = 1 or exp(F), exp{K) 1) arc 

Aj=i±\ = i — - — (— i)"' / dr exp(— mr) (mr + 1) 



2j + l 



-6jj[krjmr^^ 



B^iJij-i = * {~^Y I dr exp{—mr){mr + 1) 



,27re 



2 



x{[(-l + -^)jAkr) + -hjrikr)krdp±^^^i±^ 
.7 + 1 .7 + 1 y 2j + 1 r 

+ [(1 + fjfikr) + -j'Akr)kr]^ ,J^^ '^'f }■ (3.113) 

The connection between the wave functions u^j^i-^ij and u^-j^y^-^^ - (see Eqs. (3.64) and (3.65) appears complicated, 
but specializing as in the singlet case, we find that the terms beyond the first include higher order a terms from the 
various potential. 

For the nonrelativistic wave functions we have 

^^^==i?0±i)i.M=r^"±^X(,±i)y(r). 
We also need the small component wave function 



VjQj_ ^ exp(3g) ^^ (j - 1) - 2Qm ^^^^^ 11)^] / ^ "(j-i)ij 

r m r dr V 2j + 1 r 

U + 2) + 2q„ ^ . /ITT^ 

r dr y 2j + 1 r 



(3.114) 



For the ^Po state we have 

^'^TTP^ r°° ■71+ H~ 

Ao = / drexp(-mO{-3mno(fcr)-220 +(mr + l)(jo(M + 2j'(fcr)fcr)^}, 

3 Jo r r 

^ i2ne_ j drexp(-mr)(mr + l)(2j2(fcr) +j^(fcr)fcr)^. (3.115) 
3 Jo 

and the decay rate (3.103) involves the amplitude combination 
T.p, = (Ao + 2i?(^)) 

- /"OO 

TTl / 

= — — / arexp(— TOr){(mr + 1) 
3 Jo 



X [jo(fcr) + 2j' {kr)kr + 4j2(fcr) + 2j^(fcr)fcr] """^"'^ - 3jo{kr)mr^^^2^} , (3.116) 

r r 



in which (from Appendix D we find for this state that $++ cancels with the remaining portions of J') 

"no _ "no 
r E r 



(3.117) 
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We also have 



and 



^000 _ cxp(g) 2 d ,u+o 



E r dr 



(3.118) 



^^110Wlr=0 = ^''XlloW|r=0 = Xllo(O)- (3.119) 



The NN condition (3.41) is 



oo 



2 ^_s2 ^ + ^2 
,2r/ "-110 \ I / ""no \ I M'OOO 



ZiiM + ™ + pi<>u ] = i. (3.120) 



Our multicomponent results uses these relations in Eq. (3.116). In Appendix F we present the details that allows 
us to obtain the result of 

r(^Po-27) = ^, (3.121) 

We point out there that in the limit in which the variation of the positronium wave function is neglected (the 
nonrelativistic approximation and single component result) we obtain vanishing amplitude in the ^Pq case. As 
stressed in [40] the inclusion of the small components of the wave functions is essential for this decay. 

For the amplitude j ^ 2, I = j — 1 = 1, j' = 0, 2, 4. The relevant amplitudes arc (ignoring angular momentum 
coupling) 

27rp2 r°° [9ir (r') ?;+ (r\ 

A2 = -i^- I drexp{-mr){{mr + l){-j2{kr)+f2{kr)kr)\l- ' - ij2{kr)mr (3.122) 



and 



2jj-g2 r'x> 

= —i / c?r exp(— mr)(mr + 1) 

3 Jo 

M[lMkr) + lj'„{kr)kr]^^^^}, 



2 r°° 



b'^^ = +i — — I dr exp(— 7nr)(7nr + 1) 



x{[^j2(fcr) + i,^(fcr)fcr]^^^}, 



B^'^^ ~ —i — — j dr cxp{^mr){mr + 1) 



2 /-oo 



Mil-Uikr) + ljiikr)kr]^^'-^}, (3.123) 



with the neglect of orbital mixing where 



i,]^]^2 / ""112 12 / / mr 



mr + a r 5m"^ V V + a 



3 



1 , mr + a , 1, I mr mr + a. d.wTi-,- 



V rnr r V mr + a V rnr dr r 



W202 1 // rnr \ , V ™^ mr + a d, /2u2i; 



TO V V mr + a J r V 'Tir dr \ 5 r 



(3.124) 



Using the above expressions for V2Q2 and Mii2(?') with 



-'202 

""112 



i?ii2(r) =rxn2(r), (3.125) 
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and Eq. (3.104) leads to (see Appendix F) 

r(^^'2-27) = ^, (3.126) 

and the ratio r(3Po -> 2-i)/T{^P2 -> 27) = ^. 

Even though our approach leads to the earlier results it is of interest to see how our constraint formalism based 
approach differs from other approaches. We first note that in the constraint approach, the general frame form of the 
CM amphtudc of Eqs. (3.6) and (3.5) is 

j d^pTrT{p_,p+;h,k2)Sip ■ P)^ip), (3.127) 



in which 



r(,_,p,;fc,,fe.)=e^[Te("^) ?-";t:^^ -e("-)+7-e^- C"";^";:'l^ -^^"-^]- (3-128) 



(p_ — fci)^ + (p_ — k2)^ + m? 

In the constraint approach, from Eqs. (2.6) and (2.13) 

P 



P+ = ^-P- (3.129) 
The CM form is seen to follow directly from this since there we have p = (0, p) and 

p- = (^>p) ;p+ = (^'"P)' 

fcl = (|,k) ;fc2 = (|,-k). (3.130) 

This interpretation of the amplitude follows directly from the constraint formalism and is distinct from that used in 
other approaches which assume an on shell form for the amplitude (see e.g. [16] which uses p° = ^/m^^\^pf). The 
amplitude wc use incorporates an off-mass-shell assumption which is true for constituent particles of the bound state. 
The constraint modification of the off-mass-shell amplitude in addition places it on energy shell. This gives us the 
Yukawa modification seen in Eq. (3.12) not appearing in other approaches. 

2. Meson Decays 

a. 77^, 77^/ Decays For the j]^ the state vector is 

h,) = ^5 E 1^^)' (3-131) 

with the charge of the charmed quark equal to 2e/3. Since the interaction is color independent the resultant amplitude 
is 



9 



4V3e2 



27ri / drexp(-mr)(ji(kr)eyii>(F - K)-^](l + mr) 
^ Jo r 

v+ 

+mr^^{ji{kr) exp{F)[cxp{K) + 2sinh(A')] 
r 

+krexp{F + K)[j[{kr) - Ukr)]}) . (3.132) 

In (3.132) we take numerical wave functions from the work of [7]. The remaining parts of our multicomponcnt wave 
functions are 

^ = — ^ (3.133) 

r E r 
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and 



^110 _ cxp(g) d L' Uqqq 



J, (;r + T)— ' (3.134) 

r L ar 2 r 

which appear in that equation satisfy the TBDN condition Eq.(3.54). In the spin singlet state this is (see Appendix 
C) 



2 



/•CXD 

/ <irr2exp(2F)(exp(-2/v)[ 
Jo 

_2«;2|^-exp(-2ir)[^"»° 

,,2 o)v^^o / "000 \ I / "000 



+4«;2^ exp(-27^)[2 ( ^ ) + ( ^ ) ]) = 1. (3.135) 



The NN norm condition is 



d^r^[pM> + pioo + riio ] = i, (3.136) 



In that same hmit our NA (3.132) becomes 



9 Jo r r 

The multicomponcnt forms given by the TBDA and TBDN in (3.132) and (3.135) respectively give a decay rate of 
9.18 keV, while that obtained from the corresponding NA and NN forms (3.137), and (3.136) is 9.15 keV. If we further 
ignore the small components in these latter forms by taking Uqqq = "^oao ^-nd u^j^q = 0, then the decay rate is 9.09 keV. 
These are to be compared with the observed rate of 7.44 ±1.0 keV. Including first order QCD radiative corrections 
[43] damps these decay rates by a factor of (1 + as/n{Tr^/3 — 20/3)) giving us 6.20 and 6.18 keV.(and 6.14 keV when 
ignoring small components). For the 77^ our results are 4.81 and 2.79 keV (and 2.68 keV) respectively compared with 
the observed rate of 1.3±.6 keV. The QCD radiative corrections reduce these to 3.36 and 1.95 keV (and 1.87 keV). 
The overall additional effects of using the TBDN and TBDA above that of the NN and NA appear to be very small 
for the but for the rj'^ they arc substantial ( but in the wrong direction!). It is of interest to trace the origin of 
these contrasting behaviors. The square root of the norm (starting with a normed UQQQ/r) for the rj^ in the TBDN 
and TBDN case is 1.64, compared with 1.03 in the NN and NA case. The respective raw decay amplitudes (with the 
norm effects taken out) are 0.252 and 0.160. These are both substantial differences. However, including the norm 
effect in the amplitude cancels out these differences giving us about a 0.155 amplitude in both cases. This cancelation 
hides the substantial effects of both the TBDA and TBDN. Things are different in the case of the 77^. There the 
square root of the norm in the TBDN and TBDA case is 1.22, compared with 1.008 in the NN and NA case. The 
respective raw decay amplitudes (with the norm effects taken out) are -0.138 and -0.087. Unlike the case of the rj^ the 
effect of including the norm in the amplitude does not cancel out these differences giving us about a 0.113 amplitude 
in the first case and a -0.086 amplitude in the second. The ratio of the TBDA to the NA are 1.58 in both cases. 
However, the square root norm ratios are quite different, being 1.59 in the case of the i]^ but only 1.21 in the case 
of the 77J,. This may point to a limitation of the linear confining model used in working out the wave functions near 
threshold for the rj^, decay. 

b. Xo Decay The '^Pq decay amplitudes are from Eqs.(3.92) and in the combination from (3.103) 

^xo = ^(^0 + 25^) 



4^/3e2 27ri 



poo 

/ dr exp(— mr) exp(i^){(mr -1- 1) 
Jo 

"noW 



x{cxpK[jo{kr) + 2j'g{kr)kr + 4:j2{kr) + 2j'2{kr)kr] + 2jo{kr) sinhif}- 
-3exp(-A-)jo(M™^^^^^}, (3.138) 
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with the same color and flavor factors as before, in which (see Appendix D) 

(3.139) 



^110 _ ^110 



r E r 
We also have 

^ = ^[2-^^' + -^]<o. (3.140) 
r E r 2 dr r 

Our multicomponent TBDA results uses these relations in Eq.(3.138). The TBDN condition (3.62) becomes (see 
Appendix C) 

] r / /+\^/-\^/+^^ 



J drr^ cxp(2F) ^ exp(-2A'){ 



2./ \ \ r J \ r ) \ r 



+{2i 



110 1 I "-no 1 _ f "000 
r 



while the NN condition of Eq. (3.41) is 



+8z^2^sinh2Kf^ ) 1=1 (3.141) 



,^ ,^.2 ^_^2 ^^^2 

1 / . Or/ U 



In that same limit our decay amplitude (3.138) becomes the NA 

■^Xo ~ — n o~ / d?'exp(— mr){(7nr + l)[7o(fc7') + 2jQ(fcr)fcr + 4j2(fcr) + 2j2(fcr)fcr 



4^g2 27ri /"^ ii~ fr) 

/ J ^ -r){(77tr + l)f7o(fcr) + 2jo(fcr)fcr + 4j2(fcr) + 2jo(fcr)fcr] ^ 

r 

3exp(-i^)jo(fcr)mr^-22^}. (3.143) 







Our multicomponent TBDA and TBDN result from (3.138) and (3.141) is 3.90 keV, while that obtained from the 
corresponding multicomponent NA and NN result from (3.143) and (3.142) is 3.28 keV. If we further ignore the small 
components in these latter forms by taking u^^q = mJi^q and Vqqq = 0, then the decay rate is 0.646 keV. These are to 
be compared with the observed rate of 2.6 ±0.65 keV. The QCD radiative corrections [43] modify these by a factor 
of (1 + Q;s/7r(7r^/3 — 28/9)) to 3.96 and 3.34 keV (and 0.656 keV). The multicomponent effects are substantial even if 
we do not include the effects of the TBDA and TBDN. Those additional effects are small compared with the effects of 
including the multicomponents by themselves. This parallels that which occurs in the ^Pq positronium decay where 
the amplitude vanishes without the multicomponent (small) parts of the wave function. 
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X2 Decay The "^^2 decay amplitudes (3.92) appear from Eq.(3.104) in the separate combination 



4^/3 



and 



in which 



X2 



[A2 



2B. 
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1^202 



2(mr + l)sinh_ft:[- 
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312 
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'"ii^].72(fcr; 
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2Mi12W^ 



2 / I 
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with $ i.,A„irn, Gmp given in Appendix D and 
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with the expression for $h , App^ .., Gpm also in the Appendix D. The other radial wave functions are 
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Only the first of these latter two wave functions contributes to the decay amplitude. All wave functions contribute to 
the TBDN condition (3.62) which has the form (see Appendix C) 
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Our multicomponent results uses these relations in Eq. (3. 144, 3. 145). The NN is 
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In that same limit our decay amplitudes (3.144,3.145) become the NAs 
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Our strong potential, multicomponent result from (3.144), (3. 145), and (3.150) is 1.43 keV, while that obtained from 
the corresponding weak potential forms of (3.152) and (3.153) is 0.836 keV. Our multicomponent TBDA and TBDN 
result from (3.144), (3. 145), and (3.150) is 1.43 keV, while that obtained from the corresponding multicomponent 
NA and NN result of (3.152), (3.153), and (3.151) is 0.836 keV. If we further ignore the small and tensor coupled 
components in these latter forms by taking U112 = '^^112 ^'^'^ ""■312 '"312 ~ ''^202 = ^212 ~ 0,then the decay rate is 
0.033 kcV. These are to be compared with the observed rate of 0.528 ±.09 keV.The QCD radiative corrections [43] 
modify these by a factor of (1 — 16Q!is/7r) to 0.743 and 0.435 keV (0.017 kcV). Full tensor couplings are included in 
the first two results. As with the ^Po decay the NA and NN multicomponent effects are substantial even if we do not 
include those of the TBDA and TBDN. Those effects are themselves significantly larger than the effects of the NA 
and NN. 

d. 7r° Decay The 7r° state vector is 



1 



c—r,g.b 

where the charge of the u is +2e/3 that of the d is — e/3. Thus, the amplitude for its annihilation is modified by a 
factor of \/3[(2/3)^ — (— 1/3)^]/\/2. Otherwise the wave function discussion is the same as in the section on 77^ decay. 
So we obtain 
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Otherwise the norm and amplitude discussion is the same as in the section on 77^ decay. Our multicomponent TBDA 
and TBDN result from this is 24.7 eV, while that obtained from the multicomponent NA and NN result is 94.4 eV. 
If we further ignore the small components in the weak potential form by taking itQQQ = "i-i-ooo ^'^^ "'^no ~ ^-i then the 
decay rate is 89.5 eV. These are to be compared with the observed rate of 7.72±.04 eV. QCD radiative corrections 
modify these to 8.73 eV and 33.5 eV (31.5 eV). The influence of including the TBDA and TBDN multicomponent 
effects in the norm and the amplitude are substantial when compared to that of including just the NA and NN effects 
and bring our pion decay rate reasonably close to the observed rate. 
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e. TT2 Decay For this spin singlet decay the relevant amplitude is 
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together with the normalization condition (3.54) (see Appendix C) 
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Our multicomponent TBDA and TBDN result from this is results is 180 eV while the NA and NN result is 142 eV. If 



we further ignore the small components in the weak potential form by taking u 



202 



-''202 



and 



'112 



'112 



0, then 



the decay rate is 4.66 eV. The experimental situation is unclear. Earlier results had very large widths on the order 
of 1 keV. In the latest compilation, one result is listed as <70 eV and one at <190 eV both at the 90% confidence 
level. In any event, the multicomponent effects here are substantial which ever result we use. The difference between 
the results are small compared with the effects of including the multicomponents by themselves. 

/. "'^£'2(3872) Decay The quark-content of this state is unsure. If we assume this is a spin singlet, then 
the relevant decay amplitude has the same wave function structures as with the tt2 except for the flavor factor. Our 
multicomponent TBDA and TBDN result from this is results is 65.7 eV while the NA and NN result is 73.2 eV. If we 



further ignore the small components in the latter form by taking u 



202 



^202 



and 



'112 



ai2 



= 0, then the decay rate 



is 168. keV. The experimental situation is unclear. Using ratios given in the latest table (for the (cc, I = 2, j = 2 state) 
we take the observed value to be 435 eV. As with the it2 decay, both multicomponent effects here are substantial. 
The difference between them are small compared with the effects of including the multicomponents by themselves. 
However, unlike the 7r2 the effects are to reduce rather than enhance the rate. 
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g. 02 Decay Except for the quark content (same as with ttq, 772) this particle has a '^P2 decay amphtude and wave 
functions given as with the X2- muhicomponent TBDA and TBDN decay rate result is 31.5 keV reducing to 9.02 
keV when QCD radiative effects are included. The corresponding NA and NN rate result is 10.9 keV reducing to 
3.12 kcV when QCD radiative effects are included. The observed rate of 1.00 ±.06 keV. Including TBDN and TBDA 
effects in the norm and the amplitude are substantial and produce too large a decay rate. 

h. /2 Decay With an ss quark content, this particle has a '^P2 decay amplitude otherwise similar to that of the 
above 02. Our strong potential, multicomponent decay rate is 2.36 keV reducing to 760 eV when QCD effects are 
included . The corresponding weak component rates is 1.08 keV reducing to 348 eV when QCD radiative effects are 
included. The observed rate is 81 ±9.6 cV. As with the a2, including TBDN and TBDA effects in the norm and the 
amplitude arc substantial and produce too large a decay rate. 



IV. DISCUSSION AND EARLIER RESULTS 



A. Charmonium 



The table below (units are in keV) compares our results (both the ones that come from TBDA and NA multi- 
component results) with a variety of other quark models (ones that have not yet been subjected to the tests imposed 
on the TBDE and which for the most part do not include the light mesons in their spectroscopic calculations). 

Expt TBDE-TBDA TBDE-NA [16] [44] [45] [48] 

77^(liS'o - 2976) 7.4±1.0 6.20 6.18 4.8 10.94,10.81 5.5 3.50 

77^(2iS'o - 3263) 1.3±0.6 3.36 1.95 3.7 - 1.8 1.38 

Xo(l^Po -3415) 2.6±.65 3.96 3.34 - 6.38,8.13 2.9 1.39 

X2(l^f2 - 3556) 0.53±. 09 0.743 0.435 - 0.57,1.14 0.50.440 

Ackleh and Barnes [16] independently developed a similar approach to the one we developed for positronium decay 
[15] and then applied it to spin singlet quarkonium decay into two gammas. As in our approach, they include the effects 
of the bound state wave function on the initial decaying particle. Gupta, Johnson, and Repko [44] follow a similar 
approach. The two numbers displayed in their column correspond to two distinct approaches used in incorporating 
off shell effects. The first is similar to that used in [16] where the energy factors which arise from the Feynman 
amplitude are treated on mass shell [E ~ + m?] energy conservation, which would have E = w /2/\s not used) 
while the second set of numbers come from treating the particle on energy shell {E = 111/2 but with m? = E'^ — P^)- 
Our approach is different from both of these in that it is on energy shell, E = 111/2, but with m? 7^ E'^ — p^ it is 
off-mass-shell. The energy factors that appear in our equations are those required from the way in which the constraint 
formalism eliminates the CM relative energy- see Eq. (2.12, 2.13)). (See also our discussion in our section on ^Pq 
positronium.) To be more explicit, the portion of the Feynman propagator (p_ — fci)^ + — iO in the approaches of 
[16] and the first of [44] is treated as 2p • k±w-\/p^ + m^, in the second of [44] as 2p • k+w^/2 and in the constraint 
approach as (p — k) -\-im? . The treatment of the spin-dependent aspects of the wave function in [16], [44], and [45] is 
similar to that appearing in our earlier paper on positronium decay in [15]. We point out, however, that in our paper 
here, the spin dependent aspects of the wave function do not arise from the free spinor factors in the Feynman decay 
amplitude, but rather from the multicomponent structure of the interacting TBDE. The treatment appearing in [45] 
uses a quasipotential wave equation that gives a Schrodinger-like equation for the bound states. Their amplitude 
treatment is otherwise similar to that of [16] except that they include (as we do here) QCD radiative corrections. 
Another treatment is that of [46]. They list decay rates of 5.5 and 2.1 keV for the r}^ and r}'^ respectively, similar to 
the results of [45] (see also recent result of [? ] ) . Their treatment of the spin-dependent aspects of the wave function 
appearing is more like ours except that they use the Salpeter truncation of the Bethe-Salpeter equation but with 
energy factors in the amplitude treated on shell as in [16]. The treatment in [48] is similar to that of [46] except 
that it involves the four-dimensional Bethe-Salpeter amplitude constructed from the Salpeter solution. As with [46] 
they use a combination of scalar and time-like confining potentials. Unlike [46] and [45], and like the treatment of 
[16] and the present one, [48] then goes on to treat the light quark pseudoscalar decays. For the ^152(3872) particle, 
Ackleh and Barnes obtain a result of 20 eV close to ours of 27 eV. None of the other authors include this particle. 



B. Light Quark Mesons 



Our formalism at this stage does not include the effects of flavor mixing and consequently we do not compute the 
rates for 77, r/ 27. This leaves us with the 2 7 decays of 7r°, 712, 02, /2- We present the results of the decay width 
in the table below, including those approaches above that give predictions for some of these decays. The units are in 
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eV. 

Expt. TBDE-TBDA TBDE-NA [16] [48] 
TToOS'o - 0.135) 7.22±.04 8.73 33.5. 3.4^6.4 3.81,5.07 

7r20-D2 - 1.670) <70,190 181 142 110^270 73.2,129 

a2(3P2 - 1.318) 1000±60 9020 3120 - 766,900 

/^(3P2 - 1.525) 81±10 760 348 

Our pion rate is comparable to the others. However, the assumptions of [16] are quite different from ours. First they 
use a non-relativistic potential model for the wave function. As pointed out by [48] , [49] , and [50] standard approaches 
to the pion decay fail miserably for such models, typically too large by three orders of magnitude (by comparison our 
result is only off by 20%). Ackleh and Barnes, however, included, as did Hayne and Isgur [51] in an earlier paper, 
a phenomenological resonance mass factor motivated by an effective field theory Lagrangian \g4>F^^F^^) which 
greatly suppresses the "bare" rate. The approach we have taken above did not include such a factor. Of course, as 
they point out, the factor implied by that effective field theory is not contained within the positronium-like model that 
they and we use. (The range of values in their column correspond to a range of assumed constituent masses). The 
approach taken by Miinz [48] is much closer in spirit to the one we employ. He uses the framework of the Salpeter 
equation for the formulation of two photon decays and finds that including rclativistic effects, and the negative energy 
components of the wave function, is important even for heavy quarkonia. In addition, unlike our approach, which in 
the CM frame would have momentum space wave function dependence only on the relative three momentum, he works 
out a decay matrix element which includes relative four-momentum dependence (including relative energy dependence 
in the CM) . It is his claim that in this way, not only does the amplitude depend on off-mass-shell annihilating quark 
pairs (through the wave function) but also the exchanged quark within the diagram that arc both off mass shell and 
off energy shell. In contrast to our Eq. (3.6) his amplitude involves an additional integral over the relative energy. 
However, he finds it necessary to introduce a cutoff factor in his spectral analysis for the one-giuon exchange. In 
addition he finds that he must assume not only a different confinement mechanism for the light and heavy quarks, but 
different confinement strengths. The spectral results we obtain do not treat the heavy and light quark bound states 
differently. Further, in the two models that he considers, he finds that it is not possible to formulate the one gluon 
exchange gauge invariantly and so uses the Feynman gauge in one parameter set for a semirelativistic model and the 
Coulomb gauge in the other. By contrast, the constraint approach displays gauge invariance, and, for simplicity uses 
the Feynman gauge. 

There are other approaches that develop formalisms with natural suppression of the tt" 27 width in the quark 
model. Guisasu and Koniuk [50], using a multi-pair structure in the context of the Bethe-Salpeter formalism, show 
how the extremely bound highly rclativistic nature of the pion suppress the decay rate. The authors of [49] also 
show how the assumption of a completely diagonalized QCD Hamiltonian (with meson eigenstates predominantly qq)^ 
implies bound state effects can greatly suppress the width. 



V. CONCLUDING REMARKS 



The Two-Body Dirac equations are based on Dirac's constraint formalism and a minimal interaction structure for 
the effective particle of relative motion (first used by Todorov) confirmed by both classical [28] and quantum field 
theory [19]. This formalism displays spectral results with flavor independent interactions in very good experimental 
agreement for most of the meson spectra. At the same time, and we have stressed the importance of this in a recent 
publication [7], the formalism when treated in a nonperturbative manner naturally accounts for the perturbative 
results of QED bound states. So far this has not been fully replicated in any other approach. In a natural way it 
leads not only to good singlet-triplet ground state splittings for the light meson, but also a Goldstone behavior for the 
pion. This we showed is tied to the same rclativistic structures that account for the nonperturbative positronium and 
muonium results. Based on this and the successful off shell treatment of positronium two photon decay we had reason 
to anticipate that the quarkonium decays to two photons would be reasonable. We have found this to be particularly 
true for the ttq and 77^. There we found rclativistic effects, including most importantly the full multicomponent wave 
function and the influence of the TBDN and TBDA, to be of crucial importance. The results compare favorably 
with models based on two-body formalisms not tested as extensively as that of the Two-Body Dirac equations. Our 
pion results are unlike some of the competing approaches in that no additional effective field theory assumptions were 
made that go beyond the rclativistic potential model approach, and spectral results for all mesons are obtained in a 
flavor independent way. In light of this our results are not too unreasonable. Still one may speculate on assumptions 
made in the constraint approach which may be relaxed. For example, it may very well be that even though spectral 
results are independent of the method by which the relative time is controlled in the constraint formalism [26], the 
decay and other amplitudes may depend on this effect. The work in [26] (see Appendix A in that paper), allows one 
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to show how the relative energy restriction Eq. (2.12) which in quantum form is P ■ pip ~ or 'ip{p) = 5{p ■ P)ip(j)±) 
and this could be replaced in the amplitude Eq. (3.127) by a more general form, say 

iPip) = d{p ■ PHipi_) ^ A{p ■ P)iP{pi_) 

in which A is a distribution with non zero support [53] . In future work, having shown that the meson wave functions 
of [7] used in this paper give in most circumstances reasonable results, one will now consider applications of them to 
the meson-meson scattering process such as discussed in the beginning of this paper. 
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APPENDIX A: DIRAC MATRICES FOR THE TWO-BODY DIRAC EQUATIONS 
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APPENDIX B: RELATIONS BETWEEN FOUR COMPONENT DIRAC SUBSPINORS 

We rewrite the TBDE Eqs. (3.16) in terms of the subspinors 

0± = "01 ± V'4> 
X± = V'2±V'3- 



With Appendix A these equations then lead to 



1++ 



Eix+ - Mix. 



(Bl) 



(B2) 



37 



and so 



1 

1 



and 



X- = ~-{E2D++ + E,D++)c^^, (B3) 

Di'X- - A/i0+-i?i</'^, (B4) 

with the L'f =^ given in Eqs. (3.23a) and(3.23b). These lead directly to Eqs. (3.20) and (3.21). From Eq. (B4) we 
find using Eqs. (3.20) and (3.21) leads to Eq. (3.22). 
Using Eqs. (3.20), (3.21) and (3.30), (3.33) wc obtain 

X+ = ^^{M2[^i.p-^rT2-V(J-L + C?rTi.^2)] 
-Ml [0-2 • P - ^0-1 • V(J -L + GCTI- (T2)]}<l)+ 

- £^^{(M2 + Afi)[( p + '-V{J -L + g))- 4>^1 

+ (P + JV(J -L- Sg)) ■ + (A/2 - A/i)[ip + ^V{J - L + g)] x (f,^ ■ a 

= X+o^ + X+-^- (B5) 

Similarly 

X- = -^^^{E2[<T,-p^'-cT2-ViJ-L + ga,-a2)] 
+Ei[(T2 • p - ^0-1 • V(J -L + ga-i- (T2)]}<l)+ 

- -5^{(i?2-i?l)[(p+^V(J-i + g)).cA+l 

+ (P + ^V(J - L - 3^)) <^+o • + (^^2 + Ei)[tp + iv(J - L + g)] X 0+ • a} 
= X-ol + X--'^- (B6) 
For equal mass and electromagnetic-like interactions (J = -~g) we find 

3e;'-L' = -2F', 

-(e+i)' = -2A", (B7) 

and 

X+ - ^^{(P - ^VL) . 0+1 + [p - 1V(L + 4^;)] . ^0+0]}, (B8) 
so that with the use of Eq. (3.32) 

Xo+ = exp(F-AOr;+o = H^^(p-^VL).0+ 

= ^^P(g+^^ + -^) [p _ 1 v(L + 2^^ + 2A0] • [1 + Q„ff 

X+ = exp F(exp if + [exp(— if) — exp(if )]f f •)r7+ 

= ^{[p4v(L + 4g)]0,o} 

= 'J^2il^J£±^{[p 1V(L + Ag+2F 2if)]^+„}, (B9) 
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and so inverting the exponentials and dyads we find Eq. (3.37). Similarly, in this case, 



so that 



Xo- = 

X = exp(F + /i:)(l - (1 -cxp(-2/v))ff) -ry, 

= „ "-^P(g+^ + -^) [,p _ 1 v(L ^2F~ 2K)] X (1 + Q„ff) • ^f,^, (BU) 

giving us Eq. (3.38) again by inverting the exponentials and dyads. 

Next we work with Eq. (3.22). We transform this equation involving relations between spinors into one involving 
relations between the scalar and vector wave functions tJj^q, and tJj^q, 'tp+ ■ We begin by transforming the portion 

- E^D^+)x+ 
= eyi^{g){E2[cTi • p + ^0-2 • V(- J -L-Gai- 0-3)] 

1 

-El[(T2 ■ P + -CTi • V(- J - L - QCTI ■ CT2)]}X+ 

cxp(g){(i?2 + i?i)[(p + '-V{J + L + g))- x+1 

+ (P + ^V(J + L- -ig)) x+o-t] + {E2 - E,)[ip + iv(J + L + g)]x x+-<t}, (B12) 

and then next 

= exp(^;){M2[cri . p + lcr2 • V( J -L-gtTi- 0-2)] 

+A/i[(T2 • p + -o-i • V(J - L - t/o-i • o-2)]}x_ 
^ exp(e){(Af2-Mi)[(p+^V(-J + L + g)).X-l 

+ (P + ^V(-J + L-3g)) x+o-'^] + (^^2 + A/i)[ip+iv(-J + L + g)] xx„-t}. (B13) 

Thus Eq. (3.22) becomes 

{E2E1+M2M1) ^ 

4>^ = ^ ^+ 

-^-^{{E2 + i?i)[(p + \V{J + L + g))- x+ao 

+ (P + ^ V( J + L - 3^)) x+o-*^] + iE2 - E,)[ip + iv( J + L + g)]x x+-<r 
-{M2 - Mi)[(p + '-Vi-J + L + g))- x-<yo 

+ (P+^V(-J + L-3g)) x+o-<^]-iM2 + Mi)[ip+^Vi-J + L + g)]xX--<T}. (B14) 



For Q = — J and equal mass the matrix wave fmiction (f)_ becomes, with the use of Eq. (3.34), 

<P- = ^^IeaP " K)i'+o^o + exp(J- + K)[{1 - Q^rr) ■ ■ cr} 

^[P + ^V(L -4g-2F + 2K)] ■ [p - 1V(L + 4g+2F 2K)]^l;^,ao 

_exp(F+^+2g) + ^^(^ - 2^ - " 8^))([P " ^^^^ + + 2^))] 
•[1 + Q„,ff-]V'+)] -o- 

_^exp(2|+F+^j^p + ^^(-^ + 4^+2^ + ^K)] x {[zp - lv(L ~ 2F - 2K)] 
x[l + Qmff] • i/j+l • o- 
= exp(F - K)ij _o<Jo + cxp{F + K)[{1 + Q„,ff) • • cr. 



and so 



'^-^ ^ ^ 2igA/ ~ "2M#^[P + 2^^^ - 4C; - 2^ + 2/0] • [p - -V{L + 4 g+2F 2A')}V+o 



cxp(2g) ^ • f - [p + ^V(L - 8g-2F - 2K)][p - Jv(L + 2F + 2K) 



2EM ^ ■ -^i' ' y ■ 2 ' '"^ 2 

•[1 + Qmff] • ■0 + 

-[p - ^V(-L + 4g+2F + 2K)] X {[p + ^V(L - 2F - 2K)] x [1 + Q,„ff] • -0+} 



Using 



2K-2F = 4g, 

2K + 2F = 2L-2g, 

then gives us Eq. (3.39). 

As a check on our formalism consider the case of no potential, so that for the singlet we have 

Using the results in Eqs. (3.37), (3.38), and (3.39) we obtain 

m 

V-o = 7^+0' 

■n+ = ^pV'+o: 

The total matrix wave function is 

1 TTl 1 

*|s=0 = ^^V'+o(^ogi + —croiq2 + -p • ergo). 



We compare that with the free wave function 
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Compare this with the free spinor sohitions as in the singlet case. 



)(p)u("^'(-p) = (e + m + a-p)M("^)(0)w("^)(0) iay(e + TO + a-p) 

= {e + m + gicr-p)u("i'(0)u("^^(0)^igiCT2(£ + m + qicr-p), 



in which 



Thus 



Thus, 



{0f^ay = [x^(«) 0] 



iay 
iay 



u'-'^\p)v'-'^\-p) = (e + m + fyicr-p) 



.(^i) 



= [0 x^(^)^^.] • 

[O X^^'^'^'if^v] i£ + m + qi(T-p) 



-^[w(-)(p)^;(^+)(-p) + u^^^Hp)v(^-\-p)] 

{e + m + qi(T-p){qi + iq2)(Ti{e + m + qicr-p) 



4\/2m(£ + m) 



and 



u("-)(p)w("-)(-p) 

- — -(e + m + qi<T-p){qi + iq2){cri - icr2)(e + m + qicr-p), 

4m (e + mj 



and 



y(s+)(p)^(s+)(_p) 



1 



or in general 



4TO(e + m) 
1 

4v^™(e + m) 



{e + m + qia-p){qi + iq2)i-(7i - i(T2){e + m + qicr-p), 



[e + m + qi<7 ■p){qi + iq2)4>+-(T{e + ni + qicr-p) 



= 7r;^{'7i0+-cr+— i920+-cr + — 0+-[p - iqstr xp] + — — -^^cj)+-p(T-p}. 

Consider the special case of particle and antiparticlc at rest. In that case the above reduces to 

{qi +iq2)(t>+ -cr 



This agrees in essential terms with 



2V2 



-{qi + iq2), 



2V2 



provided that in that limit 

</>+ =ip^. 

But what about the general case? One can show that letting 



em 



-{b 1 ■ — app • +cpx )</)-) 



with 



(e + m)m 

i{qQ - 93) {qim + iq2£) 

^ = 1 — \ \ ' 

[e + m)e 

(e — m) 
b = 2- ^ '-q3, 



ego + mqs 
me{e + m) ' 



the expressions (B34) and (B27) are equivalent. 



42 



APPENDIX C: INTERACTION DEPENDENT MODIFICATIONS OF THE NORM 



In order to obtain the intereraction dependent modification 



dA 



(CI) 



of Eq. (3.42) we first need the matrix connection 



-^1 



:(/?! +751752) 



1 



iPi + 751752) cxpF(cosh/'Sr + sinhifSi-fI]2-f) 



= Co 



V- 



(C2) 



between the Dirac spinor solutions of Eqs. (2.18) and those of (3.25). 

The transformation between the 16 component column vector form of the wave function that satisfies our quasipo- 
tential equation (3.25) and the one which satisfies the Two-Body Dirac equation in hyperbolic form is given in Eq.(C2). 
The corresponding 4x4 matrix form is 

exp(F) 

■0 = ^ [cosh j^(V'_|_gi + V'-*92 + V+qo+V-Qs) 



2V2 

- sinhA'S • r{i'_^_ql + i'_iq2 + r]_^_qo+r]_q3)'S ■ f] 
A:*(r) 



where 



*(r)=^^(V'+9i + i'-m + V+qo+V^qa)- 

Whereas the normalization condition (3.42) in 16 component form is 

dA . 



(C3) 



(C4) 



d3x[V'^(l +4u)2/3^/32 



= y^Mw n\ vi ^^_] ^(/3i+ 751752) (^+5^1- ftA^.f) 

} 



x{A + BSi • f S2 • f) -(/?! + 751752) 



V- 



+Aw' J d'x{[ijl rjl V^] i(/3i+ 751752) (A + i?Si.fS2- f) 



X (A + BSi.fE2-f) -(/?!+ 751752) 



v+ 

V- 

i'-. 



(C5) 



since the matrix form of 

(l+4^'/5i/3,0)^ = [l + 2u,2(PiP,|^j + (751752 

^ /C*(r) + [-2w2^p/C*(r)p + 2w2^(-^ /C*(r)75 + S/C*(r)-S)] 
= /C*(r)+[2u;2^ig2/C*(r)zg2 + 2u-2^(gi/C*(r)gi + S/C^(r).S)] 
in terms of matrix wave functions the norm condition (3.42) is 

1 = J (fxTrij'^Cij ^ j cPxTr{K.^{v))^ CK.-^{r) = j (fxTr{{K.-^{v))'^K.-^{v) 
+ {IC'f{r))^[2w^^iq2lC'f(r)iq2 + 2w;2^(gi/C*(r)gi + S/C*(r).S)]}. 
Substituting Eq.(C3) and its conjugate this norm condition becomes 

1 = ^ J d^xcxp{2F)Tr^g{[coshK{i:lqi-i'liq2+f]lqo+i]lq3) 

- sinhii'S • f (V-'+gi - ip'iiq2 + V+Qo + V-q-s)"^ ■ f] 
X [coshK{il;^qi + %p_iq2 + V+qo+V-qs) 
-sinhii'S • f(V'+(7i +'ip^iq2 + T]^qo+j]_q3)'E ■ f] 
+ [cosh A'(7/'^(?i - i'^_iq2 + V^+qo+V-qs) 

- sinh A'S • f (V^^gi - V'-ig2 + 7]+(7o+??-g3)S • f] 

dL 

x{[2u;^^— ^ig2[cosh/i:(-0+gi +?A_ig2 + V+qo+V-Qs) 
-sinh/iTS • f('0+gi +ip-iq2 + ii+qo+V-qa)'^ ■ f]«92 
+2w^^— j(qi[coshX(V'_l_9i + ■(/^_i(72 + '/+9o+??-'73) 
-sinhiCS • f('0+(7i +i^^iq2 + t/+(?o+?/_(73)S • f]qi 
+S[coshJ!r(^/>_,_(7i + ■!/'-jg2 + V+qo+V-qa) 
-sinhifS • r{tp^qi +^p_iq2 + 'i]^qo+i]_q3)l^ ■ f]-S)]}} 

= iy"rf^xexp(2F)Tr^(cosh2if(V'+?/'+ + 7?t^++??L??-) 

— sinh 2/<('(/'^(T • f'i/'_|_cr • f + ip'Lcr ■ ri]j_(7 ■ f + r/^cr • rrj^cr ■ • f?/_cr • f) 
+{[2^2 — [cosh 2ii:(V'^.V+ - VAV-- - V+V++V-V-) 

— sinh 2_ft'('(/'^cr • rtp^a ■ f — t/iL"" • rt/j_cr ■ f — 77^ cr • i'q^cr ■ f+7yl_cr • h]_a ■ f )] 
+2w^^— 2{[cosh2/^(?/'^-0+ - V'LV'- + V+'>l+-'>]-'>l-) 

— sinh 2ii'(V;^|_(T • f^/'+cr • f — ijj'LcJ' ■ rt/j_(T ■ f + 77^0- • irj^cr ■ r—rj'Lcr ■ vri_cr ■ f )] 
+ cosh^ K{ip^_^atlj^-(T + tjj^_crt]: _-cr + 77^(T7]^-(t+77Lt77_-(t) 

+ sinh^ K[(j ■ fi/j^cr • fcr((T • fV-'+c ' f )-c+cr • f T/^Lf • fcr(cr • rip_cr ■ f )-cr 
+cr • rrj\_(T ■ rcr{a ■ hj^cr ■ r)-cr+(T ■ rrj^^cr ■ rcr{(T ■ h]_(j ■ r-cr)] 

— cosh A'sinh A'['i/'^<T(cr • fi/'+f • f + tjj'Lcrcr ■ rip_a ■ f + ?/+<Tcr • rri+(T-r+rj^_(T(T ■ nycr • f)-cr 
+ (cr • f-0^)_cr • r(Tip_^_ + cr ■ f-0^_<T • rcT^_ + a ■ fyy^cr • rcrrj^ + cr ■ f?7^_cr • fcr?7_)-(T]}] 
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For the spin singlet case (3.53) 

ip+ = ip+o(Jo; V'_ = tp-o(Jo; ■n+ = v+- ^- V- = 0. (ClO) 
Substitution and performing the remaining trace using matrix identities such as 

cr(cr fr7_|- (T(T f) cr = a{2ri^ ■ Ycr Y—ri+ cr) cr 
= -(2J7+ • fcr f-T7+ cr) 
-ri\ (T{2r]+ va r-r]+ (T) = -2r]l ■ rr]+ ■ r+rjl ■ r]+ (Cll) 
leads to Eq. (3.54) in the text. From that equation and Eqs.(3.58) and (3.59) we obtain the general radial form of 



+ ^ 2 



^ drr^ exp{2F){exp{-2K)[\^^j + i^^j ]+exp(2A')[| 

_2 sinh 2A-(- 1 + 1 /^)^ 

r ]] 2j + l r ^ 2j + l' 

/ \ 2 / \ 2 



+2sinh2A7 ^ZiZii + ^ZZiZii 

r \ 2j + l r Y 2j + r ^ 

/ + \2 /_\2 
95 / „..xr„ / "7O7 \ , / 



+4^^i;^exp(-2i.)[2 + ]) ^ 1 (C12) 

For the spin triplet case (3.61) 

ip+ = '^+-(t; ^_ = ip^-tr; 77+ = 77+0^^0; V- = ri^ cr. (C13) 
and performing the trace gives Eq. (3.62) in the text. From that equation and Eqs.(3.50), (3.69), and (3.70) and 



^_ ^ !WlY,,„^ + ""JIZDILY,^. (CM) 



we obtain the general radial form of 

i j drr^ exp(2F)([exp(2A')( I 



' - 2 



+ exp{-2K) 



.2sinh2i^[(-^±ili^,/I±l + ^^S^,/^)2 + /7 + r "b-i)i. /^). 



+{[2w'^[exp{2K)i(^ 



) 



r \ r \ r \ r 



exi,{~2K) 



_2sinh2i^[(-^^(^,/I±l + ^^(^ 

r ^2j + l r ^2j + l' ^ r y 2j + 1 r y 2 j + r ^ 

+4.^|^([- exp(2i.)( ' + (^) V 2 exp(-27.) ^ 

+2sinh2Al(-^^^,/I±l + ^^(^,0^)2 + (-!!(2±1)M ./I±l + !l(i^./^)2]]}Aci5) 
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APPENDIX D: RELATIONS BETWEEN RADIAL PARTS OF SCALAR AND VECTOR WAVE 

FUNCTIONS 

1. The Singlet Wave Function 

We start with scalar part of Eq. (3.39) 

^-0 = - + 2 VX] . [p - - VL] W,o. (Dl) 

Use 

[p + ^VL] ■ [p - i VL] = p2 - i V^i + i {VLf , (D2) 

and so 

^ = {i^^±^ - ^^[p^ - '-V^L + i (VL)^]}^,, (D3) 
By converting Eq. (3.25) to matrix form and taking the scalar part one can show (see [9]) for an alternative approach) 

-p'^ - [i {VLf ~ exp(-2e) - \v^L\^. (D4) 



and so 



^ {E^ + NP) exp(2£) 2 
^-0 i oz..^ ^171^^ exp(-2e)V'+o 



2i;M 2ME 
The other wave function contributing is from Eq. (3.37). 



^+ = S^£^^^{(p-ivL) + |exp(2/f)-l|(ff.p-ivi))*„ 



in which wc have used Eq. (3.51) in the form 



2. The Triplet Wave Functions 



In this appendix we evaluate the details of the wave functions below contributing to the triplet state. Beginning 
with Eqs. (3.39), (3.37), and (3.38) we find the relations between the radial parts of the wave functions for if)_,i]^Q, rj 
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and those of For simplicity of notation we define 



and so Eqs. (3.39), (3.37), and (3.38) become 

(£;2 + M'^) 



4>- 





{L + 20) 


2 ' 


B = 


(L ~ 2Q) 


2 ' 


C = 


(i + 6g) 


2 ' 


D = 


(3L - 2Q) 


2 




exp(2e) 


2EM ^ 



(D8) 



2£;m 

[p - iVA\ X {[p - zVS] X [1 + Q,„ff] • 
+ [p - iVC\ [p - i Vi?] • [1 + Q,„f f ] • -0+ 

»7- = -^^^[1 + Qpff-][*P + VB] X [1 + Q„,ff]^+. (D9) 
a. The Wave Function 

In this section we consider ■0_ and in the subsequent sections r\j^Q and r\_. The first portion of the first term of 
'(/'_ involves 

[p - iVA\ X {[p - iVB\ X [1 + Q„ff] • -0+} 
= [p - iA!v\ X {[p - iB'r\ X [1 + Q„,ff] • 0+}. (DIO) 



The inner portion is 



The first part of this is 



and the second part is 



Thus the inner portion is 



Hence, the first term involves 



[p - iB'v\ X [1 + Q„iff] • -0+ 

p X [1 + Q,„ff] • - iS'f X [1 + Q,„ff] • 1/;+. (Dll) 

P X [1 + Qmf f] • -0+ 

= pxtA+-^Lf (D12) 

r 

-iB'rx [1 + Qnivr] ■ ip^ 
= -iB't-xiI)^. (D13) 



[p - iB'v] X [1 + Q,„f f] • 



r 



+ 



p X -0, - — Lf -0+ -iB'fx^/>_^. (D14) 



(p — i^'f) X (p X -0^ — ^^Lf • -0^ — iB'rxxp^) 
— P X (P V'+ ^ ^^Lf • ■0_|_ — iB'rxxp^) 

-iA'rx (p X -0+ - • V'+ - iB'rxiP^). (D15) 



The first part is 



p X (p X ipj^ — ^^Lf • xjj^ — iB'vxxp^) 

^ PP ■ - P-^, - + ^K'^M-2K) ^, X . - gllip X . 

— _B"f X [r'Kij)_^) — iB'p X (rxip^). 

The second is 

—iA'rx (p X i/)^ — ^^Lf • ■(/'+ — iB'rxtp^) 
= -iA'rx (p X + ^v4'i(f X L)f • ip^ - A'B'vx (f xi/>+). 
Combining, the first portion of the first term is 

(p — iA'r) X (p X ip, — ^^Lf • ijj, — iB'rxip,) 

r 

2, .Q,n i2K~Ayexp{-2K) + A' 



r 

PP • ■0+ - P'''0+ - + — ' " X Lf • f/'H 

p X Lf • - (S" + A'B')f X (f xi/j^) - iS'p X (f x'0_^) - lA'vx (p x 1/)+). 



The second portion of the first term involves 



Qpff • (^[p - iVA] X {[p - iVS] X [1 + Q„ff] • -0+} 

= Qp{ff • p^i/j^i^ - f f • pp • ■0+ - ^5^fL • Lf • ■0+ 

fzL X f • ■0 I }. 

Thus, the entire first term involves 

[1 + Qpff] • (^[p - iVA] X {[p - iVB] X [1 + Q,„ff] • 

, 2, ,Q,n . (2g-Ayexp(-2X) + AV , 
= PP • '0+ - P 0+ - [— + ]«r X Lr • 

p X Lf • 1/^+ - (S" + A'B')rx (rxip^) - iB'p x {rxip^) - iA'rx (p x 

+Qp{+rr ■ pp • '0_|_ — ff • p"''0_|_ 2'^^ ^ ' '*^+ ^ ^ ' 0+1- 

The first portion of the second term involves 

[p - iVC] [p - iVD] ■ [1 + 0„,f f] • 0+. 



The inner portion is 



The first part is 



[p - iD'r] ■ [1 + Q,„f f] • 0+ 

P • [1 + Qm? f] • 



p • [1 + QmXr] ■ -0+ 

p • -0+ + i2K' exp(-2if )f • 1/;^ + (5,„p • f f • t/;,,.. 
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while the second part is 



The inner portion is thus 



-iD'r ■ [1 + Q,„ff] • ip 
'iexp{-2K)D'r-'ip. 



[p - iD'r] ■ [1 + Q™ff] • V+ 

p-ip++ i{2K ~ D)' exp(-2X)f • V+ + QmP • ff • -0+. (D22) 



The first portion of the second term thus involves 



[p - iC"f][p • ij}^ + i{2K - D)' exp(-2X)f • -0+ + OmP • ff • 
p[p • + i{2K - D)' cxp(-2X)f • V+ + OmP • ff • -0+] 

-iC'f [p • V+ + «(2X - i:*)' cxp(-2X)f • 0+ + <9mP • ff • 0+]- (D23) 



The first part of this is 



p[p • -0+ + j(2i^ - D)' cxp(-2A0f • ip+ + QmV ■ ff ■ -0+] 

PP ■ ■0+ + [(2is: - £))" - 2K'{2K - D)'\ exp(-2X)ff • i}}^ + {2K - D)' exp(-2i^)ipf • xj)^ 
+2A''exp(-2ii:)zfp • ff • V'.,. + Qmpp • ff • V.,., (D24) 



while the second part is simply 

-iC'f [p • -0+ + i{2K - D)' exp(-2i^)f • + 0,„p • ff • -0+] 
= -C'ifp • + C'{2K - D)' exp(-2ii:)ff • - Q„,C'ir (p • f ) f • (D25) 

This combines with the first part to give. 

[p - iC'r] [p • 0+ + i{2K - D)' exp(-2JC)f • -0+ + Q„,p • f f • 
= pp • 0+ + Q„,pp • ff • 0+ + [{2K - Df - {2K - C)'[2K - D)'\{Q„, + l)ff • ^+ 

+ (2X - D)'{Q,^ + l)ipf -1/;+ + 2K'{Q„, + l)if (p • f ) f • ^/;+ - C'if (p • + (p • f ) f • 0+). (D26) 

The second portion of the second term is 

+Qpff • [p - iC"f][p • tA+ + i{2K - DYiQ.n + l)f • "0+ + QmP • ff • ip+] 

= Qpf {f • PP • ^0+ + Qmf • PP • ff • -0 + 

+ [i2K - D)" - {2K - C)'(2/i - Dy]iQ^ + l)f • 
+ (2X-D)'(g,„ + l)if •pf-0+ 

+2K'{Q^ + l)zp • ff • -0+ - C'i{p • V+ + Qm (P • f) f • -0+)}- (D27) 

Combining the two portions gives 

[1 + Opff -llp - tVC][p - tVD] ■ [1 + Q„ff] • -0+ 
= pp • -0+ + Qmpp • f f • -0+ + Qpf f • pp ■ ^0+ + QpQmrr • pp • f f • 0+ 
+ [{2K - D)" - {2K - C)'{2K - L>)']ff • 

+ {2K - DYiQrn + l)«(pf-0+ + Qpff • Pf -0+) 

+2K'ir (p • f) f • - (Qp + l)C'ir{p ■ -0+ + Q,„ (p • f ) f • -0+). (D28) 
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Now we combine the first and second terms to give 



[1 + Qpff] • [p - -V(L + 2^)] X {[p - -V(i - 20)] X [1 + Q„ ff] • V+} 



+ [p - '-V{L + 6g)][p - '-V{3L - 2^)] • [1 + Q™ff] • xf,^ 

pp • - p ■0+ - [-^ H J«r X Lr • 

X Lf • 1/^+ - {B" + A'B')rx (f xt/)+) - iB'p x (f xi/)_^) - iA'f x (p x tp^) 

+Qp{rr • pp • •)/)^ - ff • p^V+ - ■ ■'/'+ - X f • -0+} 

+PP • ■0+ + QrnPP ■ rr■^|J_^_ + Qpvr • pp • -0+ + QpQm^v pp -vv ij)^ 
+ [{2K - D)" - {2K - C)'{2K - D)'\vv ■ 0+ 
+ {2K - D)'(Qra + l)i(pf -0+ + Qpff • pf-0+) 

+2K'iY (p • f) f • V+ - {Qp + l)C"if (p • 7/>+ + Q„ (p • f) f • (D29) 



We make the following abbreviations 



^ _ .g,„ ^ (2if-^)'(Q™ + i) + ^'^ 



and so 



r 

H = {B" + A'B'), 

I = [{2K - D)" - {2K - Cy {2K - D)'], 

N = (2A'-i?)'(0,„ + l), (D30) 

[1 + Qpff] • [[p - '-V{L + 20)] X {[p - ^ V(L - 20)] X [1 + Q„ff] • 

+ [p - '-V{L + 60)][p - '-V{3L - 20)] ■ [1 + Q™ff] • V^ 
= +2pp • 0^ - p'^ip^ + 2QpVT • pp • 0^ + Q,„pp • ff • •0_^ + QpQjnrr • pp • ff • -i/j^^ - Qpff • p^0_^ 
—Gir X Lf • — ^^p X Lf • — Hrx{rxip^) — iB'p x (rxip^) — iA'f x(p x tp^) 

+Qp{-^f^fL^f • -0^ - ^riL X f • tp^} + /ff • tp^ 

+iVi(pf-0+ + Qpff • pr-xp^) + 2K'ifp • f f • V+ - (Qp + l)C"if (p • tp+ + Qm, (p • f ) f • xp+). (D31) 
To facilitate possible cancellations expand out the triple cross products and use p • f = f • p — ^ . We need 

if X Lf • '^p_^_ = irff • pf •V'.i. — i \r\ pr■'^p_^_, 

or 

ipr-ip_^_ — iff • pf H — if x Lf • -^.i., (D32) 

and 

L X f • 0+ = (-fr • p + |r| p + 2if) • (D33) 

and 

2 

ifp • ff • -0 , = f (if • p H — )f-'0,, 

r 

Pp ff -0+ = p(f • p - — )f-0+, 

r 

if(p-f)f-V'+ = f ( if • p H — ] f • ■0_|_ (D34) 
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and 

( 2 2 \ 
f f • p (p • f ) f • = f(f-pf-pH — if • p j r • ■0+- (D35) 

Use these to simplify 

+Afi(pf -^/^^ + Qpff • pf •1/'+) + 2ii"ifp ■Tv-'il)^- {Qp + l)C"if (p • V+ + Qm (p • f ) f • 1/)+) 
+QmPP ■rr-ip^ + QpQmff • pp ff • -0+ 

= ( + ,.2 + 

+ + Qp) + 2A" - (Qp + 1)Q,„C" - . pf .^^ 

Ni 3i 
-— f X Lf • 0^ - (Qp + l)C"fzp • + Q„i{r ■ pp - — p?)?-^^ 

+QpQmr(r ■ pf • p)f • 

= — ff • + Orir ■ pf -f X Lf • ^|}^ — {Qp + l)C'fip • 

3i 

+Qm(f • PP - — p)f-0+ + QpQ,nr{r ■ pf • p)f • 0+, (D36) 



where we have used 



and defined 



pf • p = f • pp — pH — f f • p, (D37) 



M _ ^ AK' -2Q^{Qp + l)C' ^ 2QpQ^ + 2Q^ ^ 



AK' + 2QpC' _ 2Qp 



Furthermore using 



we obtain 



O = N{1 + Qp) + 2K' - [Qp + l)QmC' 

I 

= N{\ + Qp) + 2K' + QpC + + . (D38) 



f x(f X'0_^) = ff • — (D39) 



[1 + Qpff] • ( [p - ^V(L + 29)] X {[p - ^ V(L - 29)] X [1 + Q„ ff] • 0+} 

+ [p - ^V(L + 6g)][p - ^V(3L - 2^)] • [1 + Q™ff] • 0+ 

3i 

= [H - P^)'0+ + 2pp • ■j/j^^ + 0,„(f • pp - — p)f + QpQ,„f (f • pf • p)f • 1/)+ 

— (G + —)ir X Lf • tp, — ^^p X Lf • , + f — — + / ) f f • , + Ofif • pr ip, — iB'p x (rxtp,) 
r r \ r J 

Q B' 

—iA'tx (p X xp^) + Qp{+2vv ■ pp • — ff • p^'0_|_ ^fL^f • 0^ fi L x f • i/'^} 

-(Qp + l)C'fip--0+. (D40) 
Of further use are 

ip X (f x-i/j^) = — ff • tpj^ — ir ■ pxjj^ + irp-xp_^_, 

ifx(px'0_|_) = — "0+ H — ff • 0^ — if • p0^ + zpf •■0^, (D41) 
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and 



So using these we have 



if X Lf • ■0+ = irrr ■ pr-ijj^ — i \r\ pf •■0^, 
iLxr-'il)_^_ = (— ifr • p + i |?'| p — 2f) • ■j/'.i.. (D42) 



[1 + Qpff] • (^[p - '-ViL + 2g)] X {[p - |v(i - 2g)] X [1 + Q,„ff] • 0+} 
+ [p - ^V(i + 6g)][p - ^V(3L - 2^)] • [1 + Q„ff] • 0^ 
= [H - p^)V'+ + 2pp • 0+ + (5,„(f • pp - — p)f-f/)+ + QpQmr{r ■ pf • p)f • 0+ 

— {rG + N){irr ■ pr ip^ — ipv-ipj^) — ^^p x Lf • + ^ £f + f f • 

— £?'( — ij}^ — ff • xp^—ir ■ p'0^+zfp-'0_|_) — A\ — ■0_|_ + -ff • if • pi/'^+ipf •■0_|_) 

+Qp{+2ff • pp • i/).,, - ff • p^V.,. - ^^Yl?r ■ ip^ 
B' 

f (— ifr • p + i k| p — 2f) • + Ofif • pr xp, 

r 

-{Qp + l)C'Hp-xP+ 

= {R- P^)V'+ + 2pp • Xjj^ + QrnY ■ ppf ■■»/'+ + QpQmir ' pf ' P)?? ' 0+ - ^P X Lf • V + 

+gp[2ff • pp • -0+ - ff • pV+ - ^fL'f • -0+] 

+5ff ■xp^ + {A + B)'iY ■ p0_|. + Tif • pf f • ipj^ + U ipf •'0_)_ + yifp • ipj^^ (D43) 
where we use the abbreviations 



r 



T = {0-N + B' -Gr), 

U = -^ + rG + N-A', 
r 

V ^ -{Qp + l)C' - 2B' . (D44) 
Taking the above expression term by term wc need (using the abbreviations Y± = Yjm±) 

L^Y+ = [(j + l)(j + 2)]Y+, 

L2y_ = (j-l)jY-, (D45) 



and so 



, 2 d (j + l){ j + 2)\ u+ 
dr"^ r dr r 
d^ 2 d 



pV+ = - ( ^ + - '-±-^>^L^ ] r±Y^ 



dr"^ r dr r J r 



Y . (D46) 
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Also needed are 



and 



d d u+ 
«r-pi/'+ = + Y+, 

dr r dr r 

1 ^, d I—, d u+ ^ , , ^ d u+ I—, r d ^ 

fL^f.^^ = _J_[^-(^j- + l)(^^-!i^_^/7(7TT)^)Y_+jXj + l)((j + l)^ (D47) 

2j + 1 r dr r r dr r 

^ - (j + 1)—)^ + ^jU + + :r)— ]Y+ 



and 



2j + 1 r dr r r dr r 

■ - , 1 rr/jO' + l) . d .U^ r—. 7T/ + d.U>.^ 

^p-^ = 27TT^[(V^+^d;)v + v/^(^(-S^-d^)v]^- 

+ [(-^^ + U + 1);^)^ + VJijTT)i-l - ;;^)^]Y,}, (D48) 
7' dr r r dr r 

-P-^. = 2m^^^^ -^^d-r^-^^^^^^ r-d-r^^^^- 

+[U + 1)(^ + ^)^-^ + v/7(JTT)(^ - ;f 

r dr r r dr r 

1 rr/ 7(7 + 1) d^ ^U^ I— -, (7 + 1) d'^ 'U+, , 

Zj + i dr^ r r^ dr^ r 

,, 7(7 + 1) , , d^ ,7/4- 7 d^ ,it_,,, 

+ (-^^V^ - (j' + 1)tt)— + V7(7 + T)(-4 + Tt)— P49 
r dr^ r dr^ r 

- , 1 rr/i(j + l) i(j + l) .rf\it_ r— — , (j + 1) (j + 1) d d2 +,^^ 

r-ppr-V' = -— {[( "'-^ 2 -X --^Tt)— + v/.7(i + 1)(-^^^^ + ^1" + T^)— 1^- 

2j + 1 r dr dr'^ r r^ r dr dr'^ r 

+[(-^ + - (. + 1):^)^ + V7(7TT)(4 - + :;^)-]y.}, (D50) 



ff • p^i/) 



, d^ 2d (7 - l)7\,7i- r-, 7r 2d (?■ + 1)(7 + 2),7i I 



2j + 1\ \ dr^ r dr r J r ' dr^ r dr 

^/f7-^nM'^2d (j±1)(j + 2)Y^/.+ ^ ry—T rf' ^2 d {j-l)j ^u^^^^ ^ 

+ 1) (^d;:! + -rf^ )]- + vJii^)[^ + - ^— ]-}Y+ ), (D51) 

f . pf . pff . = -^[i^jf^^ + v7(7TT);^^)y_ + i-ij + 1)-^^ + v7(7TT);^— )y+], (d52) 

Zj + i dr^ r dr^ r dr^ r dr^ r 

, I r -M' , 2 J(J-1)\ 

PP-^+ = TT—ri-J 



2j + I \ dr'^ r dr r^ 

~VWTT)[(^ + '-4- - (1±^>±^] + (^^!f± + (2. + i)(. + 2) j^^^_ 



^ dr^ r dr y r dr 

1 , 2 (j- + 2)(j + l)^ 7t+ 



2j + 1 ■ \ dr^ r dr / J" 



2 d jij - 1)\ 2j + 1 d 7i_ , (2j - - 1), 7i_ - 



+ V7(JTT)[ ;^ + - - + Af:^^^:^]_}Y+, (D53) 

dr^ r dr r^ / r dr r r^ r 
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Using the simplification 

Qp+Qm+ QpQrn = 0, (D55) 

and the above identities in Eq. (D43) gives 

[1 + Qpff] • (^[p - '-V{L + 20)] X {[p - ^V{L - 20)] X [1 + Q„,ff] • V'+j 
+ [p - ^V(i + 6g)][p - ^V(3L - 20)] ■ [1 + 0„ff] • V^ 

^ ^{$__ - B'eM-20)+2Vj{ITT)^+^ _^Ar^^B^^ ^ F^d ^ ^d^u^^_ 



2j + 1 \ j,z y, r dr dr r 

+{cf>_+ + v<i(jTT)[2$++ - 26^ exp(-2g) + ^ + ^ + + + Gn^r>^]}^Y^ 

r r dr dr r 

+{-^++ + 6' cxp(-2(?) + 2V7(7TT)$_+ + % + + C,,+^^ + G,M'-^Y+ 

r r dr dr r 



+ V7(7+l)[2$- - 2^2 exp(-2C?) + + ^ + C,„, + + G,^^]}—Y+ ) . (D56) 



The quasipotentials $±± come from [9] or [34] or alternatively from taking the vector part of the matrix form of 
Eq. (3.25)) 

1^2 (j + 2)(j + l) ^ ^ 1^+^^ r2 . 

— + 5 ^*++ ^*+- — = B cxp -2^ — , 

+ + + = S^exp(-2g)^, (D57) 



in which we find 



- (^2(27TT))^'^-(^+^)^'^'+?+£(^+^) 

^3 1 2g V^L + 2) (g + i^)' 

-[(cosh2K - l)(3g' - L' - -) + (g + L)'sinh2Al^^^^^, (D58) 

r (2j + l)r 

- (^2(2^)^'^ + (2(2^ J)^'^' + X-£(-^+2^^ 
+ ^2 - ^WTT))^ + ^{30 -L- ^^) 

+ [(cosh2A' - l){20' -{L- 0)' - -) + (a + L)'sinh2/^]^^^^^, (D59) 

r (2j + l)r 

3e' L ' 



^ 2(2j + l) ^ 2 2r ^M^ {2j + \y^' ^ ' 



54 



/ 3g' L' 

\ 2 2 



1 ^3^., 3(^±^_v^g-2g'ff 



r ^ 2(2j + l) 2 ^ 2r ^ A/ ^ (2j + 1) ^ ' 

The other functions are 

= j(j + l)(Op + Qm^Y + 3^') + (2j' + J - 1)^' + e'(4j' - 2j + 2), 
C™™ = ^d^L + d^GU - 1) + i(2j + 1)L'2 - 2jL'g' + (3.7 - 1)0'^ 

Gmm = + (D62) 

A„,p = 2[(2j + l)(j + 2) + (j + l)Q„,], 

S™p = [^j(Qp-Qm)+i- ^(5 + Qm)]i' + [-3(j + l)0™ + 3jgp-2j + 7]g', 

C„p = - 39'^ - 5g'2 + 2L'g', 

F„p = [2(2j + l) + 2(g„, + Op)(j + l)-2Qp], 

G„p = L', (D63) 
App = 2.7(j + l)Q„, 

Spp = -(j + l)j(Op + Qm)(§ + 3^)' + (-2/ - 2j + l)L' - {Af + 6j)g', 

Cpp = -l^d^L + {2 + j)d^g + ^i2j + l)L'^-2{j + l)L'g' + {3j+4)g'\ 
Fpp = + 1){Q^ - Qp) , 

Gpp - jL\ (D64) 
Ap^ = 2(i-l)(2.7-l)-3jQ™, 

Bpm = [-\{j + l)Qp + ^jQrr.-j-l]L' + [3jQ„.-3{j + l)Qp + 9 + 2j]g', 

Cpra = d^L-sd^g -5g'^ + 2L'g', 

Fp„, = [-2(l + 2j)-(2j + l)Qp-2jg™], 

Gpra = L', (D65) 



For the ^Pn states we need 



and 



$++ = 4t;'2 - 2g'L' + ^ + 2V'g-^ + i^^iiiM, (d66) 

4 2 r 



Bpp = 2l' — sg' , 

Cpp = -i92i + 2a2g + i(L'-4g')', 
Z 4 

J pp yj^ 

Gpp = 0. (D67) 
Note that these coefhcients are such as to cancel the effects of <&++ in Eq. (D56) exactly! 
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For the states we require 

5 5 4 5r 5 2 r 5 

+ i^[(cosh2/v - l)(2g' - (L - gy --) + {g + L)'sinh2A'], (D68) 
5r r 



and 



and 



and 



2 " y ) ( 5 +sinli2i^) - i^_^(cosh2A- —)], (D69) 



-^4 77), 








- (Qp + Qm)(3L' + 18gO + 9^' 


' + 14^', 




2 4 




^ mm 


\1C}p 12(^771, 




Gftim. 


= 3L', 




A — 


2[20 + 3Q™], 




^mp ~~ 


[(Qp - ^Q™) - \]L' + [-9Q™ + 


6gp + 3]e', 




d^L--id^g -hg'^ + 2L'g' , 




F — 


10 + 6Q™ + 4Qp, 











= --d''L + Ad^g + -L'^ -QL'g' + m'^, 



Gpp — 2i', 



Bpm = [-^Qp + Qra-^]L'+[6Qm-9Qp + l3]g', 



(D70) 



(D71) 

5 5 4 5r 

- — [(cosh2/^- 1)(3^;' -L' --) + (g + LY smh2K], (D72) 
or r 

1 /Sg' L'\ ,cosh2i^- 1 . , orn , (^+-^)'^ , or- sinh2/^,. ,^ . 

' — - y ) ( ^ smh2A) + (cosh 2 A )], (D73) 

A_pp — 12(^777, 

Spp = -(Op + Q„0(3i + 18a)'-llL'-28^', 

ld'L + Ad^g + ^. 

2 4 

_F^pp — 12((^77), — Qp^i 



Ap-fji — 6 6(^77),, 

3 n - 1^ 

2 ^ P ~^ ^m 

Cpm = d'^L-w'^g -hg'^ + 2L'g', 

Fp,n = [-10-50p-4Q™], 

Gp^ = L'. (D74) 
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b. The Wave Function rj^Q 

For this wave function (see Eq. (3.37)) we need the single term already evaluated in Eq. (D22) 

[p - iD'r] ■ [1 + g^ff] • 
= p-xp++ i{2K - DYiQrn + l)r-'ip+ + Q,„p • ff • -0+. (D75) 



Use 



and we obtain 



with 



Thus, 



2i 

p • ff • i/)^ = ( f • p - — )f-i/)+, (D76) 



[p - iD'r] ■ [1 + Q,„ff] • 

2Q,n d_ 
r dr 



i[(2K - Dy{Q^ + 1) _ ^ _ ii]f . (D77) 



. exp(g + 2g) ^^ (j-l)-2Q„, ^^^^^r9ft' nVrn ^ni ^ 

??+o = * {[ [Qni + ^)-r + (2a - D) (g,„ + 1)J^ 



E r dr y 2j + 1 r 



[ ^^'^^^^^^'" + {Qra + - {IK - DYiQ^ + I)] p±^}Y,^. (D79) 
r dr V 2j + 1 r 

c. The Wave Function r)_ 



This wave function involves the term 

[1 + Qpff] • [ip + C'r] X [1 + Q,„ff] • (D80) 
Eq. (D14) gives the inner portion of this term as 

[ip + C'r] X [1 + Q„f f] • -0+ 
= ipxxl>^-i^Lr-x()^+C'rxiP^. (D81) 

Multiplying by [1 + Qpf f ] • gives 

[1 + Qpf f] • [ip + C'r] X [1 + Q,„f f] • V+ 
i O 

= ipxijj,+ -rL ■ Ip. - i—Li- ■ xp, + C'rxxp.. (D82) 
r 7' 



With the definition 



= -^==, (D83) 



we have 



^ ^+ ^^dr r '\]2j + l r ^ dr r ^Y2j + 1 r ^ ^ 
L.-0+ = 0, 



and so 



^+ ^ r V2j + 1 r V2J + 1 ' 



[1 + Qpff-][ip + C'f] X [1 + g„ff-]i/j+ 

• n^ ^ ~ 1 jQ'n , ^,^1 / j + 1 M-(r) 
'^^^d^^^^^ ~ + ^^V2j + l r 



dr r r V 2j + 1 r ^ 



APPENDIX E: PLANE WAVE INTEGRALS 
1. States with j = I 

For general j — I states we eonsider separately integrals from Eq. (3.75) of the form 



Consider first 



We have 



J d^rexp{-ik-r)rg{r)Yj,n{^), 
y dVexp(-zk • r)/±(r)Y,„±(f2). 

J dVexp(-ik-r)f5(r)Yj>„(rj). 



f = sin 6 (cos 0x + sin 0y)+ cos 9z 

yio(0)z - J—Yn{n)i^-^y)+J — Y,^,m^+^y). 



Here x, y, z are arbitrary miit vectors fixed in space. And so 

J d^r exp(-ik • r)fg(r)Yj„j(r2) 



-(x-zy)^^ J dVexp(-jk.r)f5(r)yii(n)y,,„( fl) 
+ {x+iy)^ J dVexp(-ik-r)f5(r)yi_i(n)yj„,( fl) 



Now with 



exp(-zk.ry)=4^^ ^ {~iY' j,, {kr)Y;,^,{n)Y,,^, {n,), 
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in which we define the angles ^lk{9k, 4'k) relative to the fixed unit vectors x, y, z, we use 



J dnY*,^,Y,^„Y,^ = yJ^y^Ol; 001/0)01; mm"|/m + m"); /+j + leven 

= 0; / +J + 1 odd, (E6) 

and find 

J dVexp(-ik • r)rg{r)Yj„r {n) 



-^^^r,.„+i(n,)(jl;ml|/m+ 1) + l^^y^.,,„_i(f2,.)(jl; m - l|/m - 1)]. (E7) 
Next we examine 

I dVexp(-zk.r)/±(r)Y^,„±(f^) 
= y"rf3rcxp(-ik-r)/±(r)[a± frj™+r6±py,„]. (E8) 



Use 



and Eq. (E7) 



Likewise 



J dVcxp(-ik • r)f±{r)pYjm{i^) 
= k J d^rcxpi-ik-r)f±ir)Yj„,{n)+i J d\ cxp{-ik ■ r)rf'^ir)Y,^{n), (E9) 

J dVexp(-zk-r)/±(r)Y,„±(r!) 

I dVcxp(-ik • r)[(a±/±(r)+i6±(r/4(r) + f±{r))rYj^+rb±f±{r)kYj^] 
ATTkYj„-,i^k)i-iy I drr^jj{kr)krb±f±{r) 

POO 

J drr^ f±{r)[{a± ~ 2ih±)jj,{kr) - z6±jj, (fcr)fcr][zy,.,„(rjfe)(jl; m0|/m) (ElO) 
- ^'^~^\ ,.^+i{nk){jl; ml\j'm + 1) + i^±|^y^.,„„i(nfe)(jl; m - 1|/to - 1)]. 

y d3rexp(+ik-r)/±(r)Y,„±(r!) 
-47rkr,,„(^^fe)(-^)■'■(-)■'■ / drr^jj{kr)krb±f±{r) 







/•OO 

; / drr^f±{r)[{a± - 2ib±)j,,{kr) - tb±j'j,ikr)kr][zY,,,n{nk){jl; mO\j'm) 
Jq 

^2^-^Y,,„,+,{nk){jl;ml\j'm+l) + i^^y,.,„_i(n,)Ol; m - - 1)]. (Ell) 



59 

We choose our unit vectors to be defined in terms of the photon decay momenta and transverse polarization vectors, 

(E12) 



i^^iy) _ (±) 

V2 



With this choice, the spherical harmonic reduces to just 



Y^mi^k) = SrnO\l'^^^. (E13) 



Then Eq. (E7) becomes 

J d^rexp{-ik-r)rg{r)Yj,n{0) 



= VMSj + T) E ^^^y M ^ O'l; 001/0) / drr^j,,{kr)g{r) 

X [kOl; 00|./0)d™o - e'^'Ol; -ll|/0)5,„_i + e^+^O'l; 1 - ll|j'0)5™i], (E14) 



and Eq. (ElO) becomes 

J dVexp(-ik • r)/±(r)Y„„±(r!) 

= ^/M2]TT)kS,nn{~^y / drr\],{kr)krb±f±{r) 

Jo 

+ + 1) J2 ^^^^^^^^O'l; 001/0) / drr'f±ir)[ia±~2ib±)j,,ikr)~tb±f^,ikr)kr] 

x[k(jl;00|j'0),5,„o - e(-'0-l;-ll|/0)5„„i + eW(jl; 1 - ll|/0)5„a]. (E15) 

In the expression Eq. (3.75) the prefactor e^^^^ x e'"^) is orthogonal to e'^*^ Hence, only the k terms in Eqs. (E14) 
and (E15) will contribute and including forms for g{r) and f±{r) 



/7re e^^' x 



-exp(-i^)( "^^-^^^^-^'^ Y_ + !W,(!l y^^^^^^pHH] 
r r r 



27re2e("^) xe("^).ky(27TT)( (-i)-'" 00|/0) (jl; 00|/0)<5™o 



X drr^ exp(F + K) |^£^P(_JZ!!l^ ^ (fc^ 



j'=\j-M 
exp(-mr)V^. ."jOjW 



-S^oi-^y I rfrr^m""P("'"^^,(fcr)fcrexp(f-X)(6/^+^^^^-^'^+6_''^-^^^^-^'^ 







- E (-0^'^^^^^^(jl;00b-'0)(jl;00|j'0)<5,„o 

X drr2m^^^t!!^[(a^ _ 2ih+)jr{kr) - ib+j'{kr)kr] exp(i^ - K) "U+m^'^'> 
Jo r r 

+ [(a_ - 2ib^)jy(kr) - ib^j'y(kr)kr] exp(i^ - /f) '^"^^^^ ] j , (E16) 
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and with Eq. (3.52) this leads to 

^/^e2e("l) x e("=) • J dV[cxp(-i k • r)- cxp(+ik • r)] 



\ r 



x( (~'y'(^ (-)' V^ ' ^ jl; 001/0) 0-1; 001/0) J„,o^°°drr^ (^£^Ph!^y ^ (fc,) exp(^^ + if)>^-^'^ 

Jo r i + 1 

^ (-*)^'(jl; 001/0) (1 - (-)^')ii^l-ll!^0-l;00|j'0)<5,„o 



/o ' J- + 1 y 2j + 1 r j\l 2j + l r 



X / drr 



m 



r .7 + 1 J + 1 V^J + l 



+ [(1 + |).?,,(fcr) + i/,(fcr)fcr]^^^^cxp(^^ - j , (E17) 

and to Eq. (3.80). 

2. States with j = Z ± 1 

We seek from Eq. (3.84) the evahiation of the dyadic integral of the forms 

j d'^rexp(— «k • r)G(7')ff, 

j d3rexp(-ik- r)F±(r)Yj„±(rj)f 

= y"rf^'exp(-ik-r)F±(r)[(a±f+&±rp)Y,„(f^)]f. (E18) 

Consider the portion 

j dVexp(-ik • r)i^±(r)[pY,„ ( n)]v. (E19) 

Integration by parts gives 

j dVexp(-ik-r)rF±(r)[py,>„( f2)]f 
= j d^rYjrriifl) exp(-ik • r)[kfrF±(r) + iF±{r)l + irrrF!^{r)]. (E20) 



Its trace is 



J dVYj„(n)[k • rrF±{r) + 3iF±{r) + irF!^{r)] exp(-ik • r) 

/d 
dVr,™ (J^)[F±(r)i— + 3ii^±(r) + irF4(r)] cxp(-zj/k • r)|^^i 

An{-iyYj,n{nk) / r^dr [F±{r)t— + 3iF±{r) + irF:^{r)]j,{ykr)\ 

Jo uy 



= 4:Tr{-iyYj,n{^k) / r'^dr [F±{r)ij'j{kr)kr + 3iF±{r)jj{kr)-3iF±{r)jj{kr)-irF±{r)kj'j{kr)] 
Jo 

= 0. (E21) 
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Hence the trace 



47r(-i)^yj™(f2fc) / r^drjj{kr)a±F±{r), 



and with Eq. (E13) leads to Eq. (3.89) upon using Eq. (E13). Continuing with the dyad, note that 

J d^rYj„,{n) exp(-ik • r)krrF±{r), 



(E22) 



(E23) 



will give zero contribution to the amplitude Eq. (3.84) since it is sandwiched between transverse polarization vectors. 
The remaining part is 

/•oo poo 

4Tri{~iyYj^{nk)l r^drF±{r)jj{kr)+i ffdV exp(-ik • r)r^^4 (r)r,„, (E24) 
Jq Jo 

which because of the integration by parts is symmetric. And so, effectively we have 

J d\ exp(-ik • r)^^± {r)Yj^± {n)r 

= J dVexp(-ik-r)[a±F±(r)+i6±rF4(r)]ffY,-,„+6±47ri(-i)^Y,„,(f2fc)l^ r^drF±{r)jj{kr)]. (E25) 

Since this dyad is sandwiched between polarization vectors transverse to k we need only consider the portions of 
ff and 1 orthogonal to k. Hence we replace 1 by 1 — kk and ff by (1 — kk) • ff • (1 — kk). As in the above section 
we choose the unit vectors in Eq. (E12). So, with 

(1 — kk) • ff • (1 — kk) ~ sin6'( cos(/)x+sin(/)y) sin0( cos(/)x + sin0y) 



3 V 5 

V 15 V 15 



(E26) 



we find 



(1 - kk) - y d-'rexpi-ik ■ r)G'(r)ffy,,„ -(l - kk) 
= (e(+)e(-)+e(-)e(+)) y d'V exp(-ik • r)F±(r)i [Vi^Foo - y^>2o]>S. 
J d^rcxp{-ik ■ r)G(r)r2-2iGm 
+e^-)e^-^^ J d^r cxp{-ik ■ r)G{r)Y22Y,^. 
Using Eq. (E5), Eq. (E13) and 

J dnY;,„,^ 



(E27) 



i2j + 1)5 



y 47r(2.7' + 1) 
0; j+/ + 2 odd. 



(j2; 00|/0) (j2; mm"\j'm + m") ; j + j' + 2 even 



(E28) 
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we have 

(1 - kk)- J d^rcxp{~ik ■ r)G(r)ffY,,„ • (1 - kk) 

^ j' = \j--A 

(^^ + ^)(6(-)e(-))4^y| E ^^^^^^^(j2;00|/0>0-2;2-3|/0)<5„,2 j r^dr{-iy' j,,{kr)G{7{^2<d) 



An 
and so 



J' = |j-2| 



(l-kk)-y rfVexp(-ik-r)i^±(r)Yj„±(r2)f • (l-kk) 

(1 - kk)- y" dVcxp(-ik • r)ffYj„[a±i^±(r) + ib±rF!^{r)] ■ (1 - kk) 

+ATTib±{~iy{e^+'>e^-'> +e^-'>e^+^) / r^drF±{r)jj{kr) 



1 f 

+ e(-)e(+))-y4^(27TT){5„o / r^dri-iy j,{kr)[ia± + 3i6±)F±(r) + i6±rF4(r)] 
■J Jo 

- E ^^^^^^^^0-2; 001/0)0-2; 00|/0)<5™o H r^dri-tf j,,ikr)[a±F±ir) + ib^rF^^ir)]} 

+ E ^ii^0-2;00|/0)0-2;-22|/0)._ 

i'=li-2| 

/•oo 

X / r^dr{-iy' jji{kr)[a±F±{r) + ib±rF!i.{r] 
Jo 

+ ,.<-..<-.);E?Z±i) f " + '-"^"''' (;2.00|/0)(,2.2-2|/,)..„. 



X / r^dr{-iy jj,{kr)[a±F±{r)+ib±rF;^{r)]. (E30) 
Since these expressions involve symmetric dyads wc find together with the trace part 

-(eWe(-) + e^-)e^+)) J dVexp(-ik • r)F±(r)Yj,„±(f7) • f 
= _(e(+)e(-) +e(-)e(+))47r(-i)y^:L±l5„^oy r^drj,{kr)a±F±{r), (E31) 
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that 

-iy^e^ J d^rexp{-ik-r)F±{r)[ Yj,„±(f7) • e^^^'f • e^"^) + Yj,„±(r2) • e^^^^f • e^^i) 
-Yj„±(17).fe("i) 



2 E 2 ^ 0-2; 001/0)0-2; 00|/0)<5„o / r2drH)^'j,. (fcr)[a±i^±(r) + z6±ri^4(r)]}(e(°^).e(°=)) 



-3Smoi~^y / r^drj,(fcr)a±F±(r) 
Jo 

j'=li-2| 

/ ■ -/ 

-zW^^i^^ 5] (i±l_l)!!il(,2;00|/0)0-2;-22|/0)<5„_2 

j'=b-2| 

/•OC 

X / r^dr{-iy' jj>{kr){a±F±{r) + ib±rF;i^{r)]e^°"^ ■ (e(+'e'+)) • e^^^^ 
Jo 

-M,.'JMpl g <' + '-"'""' (,2;00|/0)(,2;2-2|/0).,.. 

j'=b-2| 



X / r2(-i)j'dr/,(fcr)(a±i^±(r) +i6±rF4(r)]e("i' • (e(-'e(-') • e("2). (E32) 
Jo 

Integration by parts gives 

-iV^e^ J dVexp(-ik- r)F±(r)[ Yj,„±(f7) • e^^^^f • e^"^) + Yj„,±(r2) • e^^^^f • e^^^) 

-Yj„±(f]).fe("i) 
.27re2 



-«^\/(27+l){(-«F(5mO / r'dr[-jj{kr)a±F±{r)-2tb±j'^{kr)krF±{r) 







-2 X] ^^i^^^^(j2;00|i'0)(j2;00|j'0)^™o 

j'=li-2| 

/•OO 

X / r^dr{-iy'jj4kr)[(a± - 3ib±)jj,{kr)F±{r) - ib±j'{kr)krF±{r)]}{€ (°i).e(°2)) 



^W^^i^^ ^ (i±i_l)!!!l(,2;00|j'0)0-2;-22|/0)<5™_, 

J' = |j-2| 



/•CXJ 

; / r2dr(-i)^'[(a± - 3ib±)jf{kr)F±{r) - ib±j'y{kr)krF±{r)]e'^°'^^ ■ (e^+^e^) • e 
.«„»^/M^ f " + '-'>"'' (,2;0O|/0)(,2;2-2|/0).„. 



3 

i'=b-2| 



/■OO 

/ r^dr{-iy' [[a± - 3i6±)jj' (fcr)i^± (r) - i6±jj, (fcr)fcri^±(r)]e("i^ • (e^-^e^"^) • e^"=\ (E33) 
Jo 
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and simplifying 

-iy/^e^ J (frexp{~ik-r)F±{r)[ Yj>„±(r2) • e^^^^f • e("=) + Yj„±(f7) • e'^^^f • e^^^^ 

-Y,„±(17).fe("i).e("^)] 
27Te^ 



z— -VW+1){(-*F / r^dr[-j,{kr)a± - 2tb±f^{kr)kr]F±{r) 







J+2 

-2 ^ (-iy'(i2;00|j'0)(i2;00|j'0)<5™o 

i'=|j-2| 

X / r2dr[(a± - 3i6±)jjv(A:r) - i&±jj,(A:r)fcr]F±(r)}(e(°^)-e(°=)) 
Jo 

i+2 



y 2(2j_+2) ^ ^2^^[(^^ _ 3i6±)j,.(fcr) - z6±j;.,(fcr)fcr]F±(r) 

[0-2;00|j'0)0-2;-22|/0),5™_2]e("^) • (e^+^eW) • e("=) 



j'=b"-2 

X 

2 

+ 0-2;00|/0)0-2;2-2|/0)<5,„2e("^) • (e^-^e^")) • e*"^)]. (E34) 



FromEq. (E29) 

-iy/^e^ J (frexpi^ik ■ r)G(r)f • e^^^'f • e^"^' 



- E 001/0) (.72; 001/0) J,„o /"r2dr(-i/j,,(A:r)G(r)](e(«^).e(°^)) 

+V6 E ^^^^^^^^0-2;00|j'0)0-2;-22|/0)5™_2 / r2dr(-*/j,. (fcr)G(r)e("i) • (e(+)e(+)) • e^"^) 

i' = li-2| 

+^ ^ il±l_lZm(j-2; 001/0)0-2; 2 -3|/0)<5™2 / r2dr(-*)^'j,. (fcr)G(r)e("i) • (e'-'e^-') • e("^)(p35) 



i' = li-2| 

Combining and including for g(r), F±{r) and G(r) 



3=l±l'^27 
2 



X 

U 



J (fr[cxp{~ik ■ r)+ cxp(ik • r)] cxp(F) ^ ^ ^^P( ™'") ^ 



. + [„exp(-A-)^ - 2(„. + l/,-)™hA-(-,/iiiS«lii + . ,(„)2St!2lj 

r \/2j + lr r\/2j + l r 
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27re2 r, 7 

-z^V(2j + l) 



({(1 + {-m-^y I r'dr { !^pL^\ exp(F + K) 



2 E (i + (-r')HF' ^^^^^^^0-2; 001/0)0-2; 00|/0)J„o 



j'=li-2| 
x[/ r'^dr 



+ 1(1 + j)fr(tr) + jj;-('--'-)*'-l|/ 2/+l 0"''"' ■ 

j'=b-2| ^0 V 7- 



exp(F + K){[{-1 + --^)j/(fcr) + -l-^f^,{kr)kr]^ 



+ [(1 + -)3r{kr) + ij;,(fcr)fcr]./-^^^^^^ili^} 

J J V + 

x[0-2;00|/0)(j2;-22|/0)(5^_2]e("^) •(£(+)£(+)) 
+ (j2; 001/0)0-2; 2 -2|/0)5™2e("^)-(e^-)e(-)) 



/exp(-mr)\ / j + 1 , "(j-i)ij / .7 



V r / y 2j + 1 r r V 2 j + 1 



E ^^^^^02; 001/0) 02; 00|/0)<5,„o(l + (-)^') 



/•OO 



V6 E ^^i^^^^O2;00|/0)O2;-22|/0)<5™_2(l + (-)^') j r^dr{^iy' j,,{kr) 



j'=b-2| 

V r J y 2 j + 1 r r y 2j + 1 

V6 E ^^^^^^^02; 001/0)02; 2 -2|/0)<5„,2(l + (-)^") j r^dr{~iy' j,,{kr) 



2 

j'=b-2| 



x[-4sinhK f £^Ph!!^) [_ /2±1!1(£±1M + /_^],(oi) . (,(+),(+)) . 

^ V ^ y ^ Y 2j + l r r Y 2j + lJ ^ ^ 

+(1 + / exp(-mr) exp(F)r/ (fcr)3[m exp(-is:)-^ 







-2(„. + Ur) .n.KiJl±L"j^ , ^J-l-,,^. . ..... 



(E37) 
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or 



{I + {~y){~iy 6rnQ I dr cxp{~mr) cxp{F) [ — 3jj{kr)[mr cxp{—K) ^^^^ ^ 







+(mr + l)exp(/i){[fe(fcr) + ^—-^f^{kr)kr)Jj^^^±l^ 



-4si„hA-(-»r+ 1)1- + J_L_|,-.,,,,),<«.. 



E 2 -{j2;0Q\j'0){l + {-y '-{~iY[2j^ drcxp(-mr)cxp(F)(mr + l) 



i+2 

E 



-exp(AO{[(-l + -^)jAkr) + -L^f^,ikr)krJ + [(1 + ^)j,,(fcr) + i,;.,(fcr)fcr] J-L^^^i^^il^} 



+jf {kr) [2 sinh K{mr + 



2j + 1 r r y 2 j + 1 

x(^(e(°i).e(«=))(j2;00|/0)5,„o 

-^/6[e("l) • (eWe(+)) • e^"^) (j2; -22|/0)(5,„_2 

~ V6e("i) . (g(-)e(-)) . ^(-2) (_^-2. 2 - 2|j'0)J™2}) (E38) 
From this we obtain the result Eq. (3.91). 

APPENDIX F: POSITRONIUM DECAY 

111 this appendix we focus on positronium decay ampHtudes assuming weak potentials (£, /C = 1). We wish to 
verify that our relativistic formalism gives the standard results for the ^Sq,'^ Pq, and decays. 

1. ^Li amplitudes 

The ^Li amplitude is 

i+i 

Mii,^27 = -v/2jTTe("^)xe("^). £{^^=,(1 - (-)^)5,„o + ^ GgCl - 00|j'0) (jl; 00|/0)<5,„o}, (Fl) 
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in which 



G'fJ = -2TTe^{^iy' J drcxp(~mr)\(mr + l)[jj,{kr) 

+ 



J J V 2.? + 1 



+ [(1 + -)jAkr) + -j'{kr)kr]J .-\ ^ ^■''^'^^ . (F2) 



We tabluate here the portions of the wave function that will contribute to the same order in a to the decay amplitude. 
We consider the relation between the various contributing components of the wave function. The relation between 
u~Qj and UjQj from Eq.(3.58) is u~g^- = ^uJqj. For electromagnetic interactions only (see [41]) 



M = m,E='-l^l-'-^ (F3) 



and for positronium A = —a/r. For positronium, w = m + 0{a^) and so 

M 



E V '^'^ + Q! 



(l + 0(a^)) (F4) 



Next consider the contribution of the small component wave functions £^nd From Eq.(3.59) for 

positronium we have 



^0-1)1.- ^ 1[ d (.7 + 1) 

r m dr r ''^^"ly 2j + I 



-[zr-z]y^.m\ i7^ + 0{a'). (F5) 



We focus on the decay amplitude for ^5*0 positronium {j = 0, j' = 1). In that case, our amplitudes are 



1 ,9 exp(— mr) ,1 d 



Fq = -2i7re / drr m ji(.kr)kr{——ip„„Q), 

Jq r m dr 

gI^^ = —2ne^ / dr cxp(—mr)((mr + l)\ji(kr)t -0 

Jo ' \ mr + a 



+mr[ji (kr) + f„ (fcr)fcr] ^^4'ooo) , (F6) 



- dr 

with the nonrclativistic wave function given by 



^000 = 7e=- exp(-amr) = (F7) 
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Because of the short ranged exp(— mr) factor, the exponential wave function can be replaced with its value at the 
origin. Hence (with k = m(l + 0(a^)) we have 

n, , (ma)^/^ f°° , 2 exp(— mr) 
Fq ~ — ziTre (—a) j= — / arr m ji{mr)mr 



Stt Jo r 



-2iTT{'iTTa){-a)- I dxx^ exp(-a;)ji(x) m^/^a''/^, 



dr exp(— mr)((mr + l)[ji(mr)w 

Stt Jo V + " 

mr[ji(mr) + j[(mr)mr][--a)) 



= — 27r(47ra) / dccexpf— a;)((a; + l)h'i(a;)w axcossh'ifa;) + (x)a;]) 

mV87r Jo \l x + a 

= -27r(47ra) / dxeM-x){x + \)h{x).r^^ + 0{m^l^a"^). (Ff 

mV87r Jo ' V ^ + 

So, the small components do not contribute to the singlet decay rate at the order m^/^a^^^. The term that does is 



^(1) (toq;)^/^ 87r^m 



g(a) = / dxji(x)exp(-a;)(l + x) J^— 
Jo \l x + a 

= .g(0) + ag'(0) + ..., 

with 

5(0) = ^ dxji(a:)exp(-a;)(l + x) = i, (FIO) 

and 

1 /"^ 1 

9'{0) = -- dxMx)cxp{~x){l + x)-, (Fll) 

2 Jo a; 

which is finite (the integrand behaves like a constant at the origin). Hence to lowest order we find 

= mi/2«5/22V8^2 ^ 1 ^ 2V2m'/^a^/W (Ff2) 



which gives the correct form for the decay rate 

G 



For later use define 



and 



Now we have 



Kj^k= / dxjj{x)exp{~x)x^^^\ 2j + fc > -1, (F14) 
Jn 



fc=0 ^ ' 
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and so 

z^i ^ ;!(2/ + 2j + l)! ^ ' 

We evaluate numerous related forms of these summations for several sets of needed indices k in the section following 
the one below. 

^ij=i±i amplitudes 
The ^Lj^i±i amplitudes are 

„2 



27re^ f°° ( 

Aj^i±\ = i (—*)"' / dr exp(— mr) I (mr + 1) 

3 Jo V 



'3jj{kr)mr- 



(■•'■) 27re^ 
BjLl_^^ = i / dr exp(— TOr)(TOr + 1) 



+ [(1 + -).,,,(fcr) + -jHkr)kr]J J^^ ^'-'^^'^ }. (F18) 
We neglect the effects of orbital angular momentum coupling. Then 



A,=i^i = ^^i~^V r drexp(-mr)((mr + l)fe(fcr) + f^{kr)kr) J p^''-^^^^^^ - 3j,(A:r)mr^^i^) , 
6 Jq J + -L y a? + ^ 

(F19) 

and 



= / drexp(-mr)((mr + l)(-j,(fcr) + -j;.(fcr)fcr)W — L^^^^i^i^ 3j,(fcr)mr^2liZ), 

3 Jo J \l 2] + i r r ' 

(F20) 

The connection between the wave functions w^.^^jj^^. and tt^-i-^jij (see Eqs. (3.64) and (3.65) )appears complicated, 
but specializing as in the singlet case, we find that the terms beyond the first include higher order a terms from the 
various potential. Simplifying we find 

_ exp(3g) ^ F„„ ^^2^) 



and 



in which 



E r m? {2j + 1) r dr 



exp(3^;) App Fpp d 



E r vn? (2j + 1) r dr 



Amm '^ji.j ^^Qrn App^ 

F^nn, = 2j{j + l){Qp - Q,„) = -Fpp, (F23) 
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where 



i = (\/zz^-i) + o(«')' 



1 — 2A/w V + ^ 



V mr 



Thus 



, 3 

mr 



mr + a r m^ (2j + 1) y ymr + a 



1 mr + a 1 / mr mr + a d '^(j-i)ij 



and 



_ / mr 2j(j + 1) 



, 3 

7Tir 



mr + a r m^ {2j + 1) V \mr + a 



x{[-l( - 1) + -(J^^- /!!!!:±i!)^]!l(z±ini}. (F26) 

r-^ V mr r y mr + a V mr dr 

For the nonrclativistic wave functions we have 



,+ 



u7_ 



r 



We also need the small component wave function 



Vjoj _ exp(35) - 1) - 2Q,„ , / J "(j-i)ij 



-{[- ^-(Q™ + i):r 



m r dr y 2j + 1 r 



For j = = Z — 1 we have 



I drexp{-mr){{mr + l){jo{kr) + 2j'Q{kr)kr)J ■ — rxiio('^) 

V + ^ 

/ f mr \ d 

3.7o(fcr)W — (Q„ + 1)[2 + r-]rxnoM) 

y \ mr + a J dr ' 

f dxexp(-a;)((a; + l){jo{x) + 2j'o{x)x)J — ■ — rxnoW 
\j x + a 



(F24) 



— - 1) + -(a/^— - V — );r]^^^}' 

mr r V mr + a V mr dr r 



27re 

= i 

3m 

-3.7o(x)^[2 + x-f]rxnoW), (F28) 
where as before we let x — mr ~ kr. As before to lowest order we can replace 1 giving us 

Ao = *^Xiio(0)^ dxexp{-x)x{{x + l){jo{x) + 2j'^{x)x) - 9jo{x)) . (F29) 

(The factor of 9jo(2;) in the integrand would be absent if we had ignored the small component v of the wave function 
ofEq. (F27).) Similarly 

= «^Xiio(0) / dxexp(-x)(x+l)x[2j2(x)+j^(a:)a;] (F30) 



or since 



we have 



where 



Integration by parts gives 



where 



and 



Since the decay rate is 



and 



and so 



and we obtain 



Since 



we obtain 



^i?110(f)|r=0 = ^rXiioW|r=0 = XlloW 



^0 = ^ 3;^d^^"°WI^-^^ 



h = j da;cxp(-x)x((x + l)(jo(a;) + 2jo(x)a;) - 9jo(a:^)) 

"'0 

poo 

jf = I (ixcxp(— + l)x[2j2(a;) + J2(x)a;] 

"'0 

/•OO 

A'j.fc / (ia;jj(a;)exp(-x)x''+''; 2j + fc > -1, 
Jo 



rC^Po - 27) = ^1^0 + 24''|2 = i^i?;^„(0)Vi + 2jP'r 



ATo.s = Aro,2 = A'o,i = - 
K2.1 = 1 



2 

Jf ) = 1 



rfPo^27) = — i?iio(0) : 



(^(n=2)llo(0) 



(24) (32) 



n^po - 27) 
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which agrees with the results of [39] and [40]. We point out that without the smaU component wave function (F27) 
we would have obtain zero for this rate. 

Next we consider the case for j — 2, — 1 — 1, j' = 0,2,4. The relevant amplitudes are (ignoring angular 

momentum coupling) 

o 2 foo PliT (r\ ?;+ (r\ 

A2 = -i—- / dr eM-rnr)i{mr + l){-n{kr) + j'(kr)kr)x - ^^"^^ ' - in{kr)mr M , (F40) 

3 ./n ■ V 5 r ■ r ' 



and 



2^g2 /"OO 

i?2°^ = —i / c?r exp(— mr)(mr + 1) 

3 Jo 

x{[^Jo(fcr) + i,^(fcr)fcr]/|^^}, 



= +i / c?r exp(— mr)(TOr + 1) 

3 Jo 

H[\j2{kr) + \Ukr)kr]sl\'^^}, 



i?2^'' = — i — - — / dr exp(— TOr)(TOr + 1) 



.27re 



x{[^j4(fcr) + i,^(fcr)fcr]./f ^^}, (F41) 



with the neglect of orbital mixing where 



12 / / mr 



mr + a r 5m^ V V mr + a 



^ 1 / mr + a -[_)+"'"( / / mr + a d uj^2 1 

V Tir r V mr + a V dr r ' 



W202 



1 // mr y mr + a d 2 uj^ 



r TO y V + a / r V "ir dr \ 5 r 

Using the above expressions for V2Q2 and Wii2('') with 

^ = Rii2{r)=rxu2{r). (F42) 

and approximations used earlier we find 

A2 = -i — - — Xii2(0) / dr exp{-mr)(^{mr + l){-j2{mr) + j'2{mr)mr) 
3 Jn 



5 V "^r + a 5to^ V V + 



/2 , I mr 12 // mr V^r 1 / hnr + a . 

xV^Ha/^^— ~— W( 3^— ) [--2(\/— T— 



1 , I mr I mr + a ^ 

mr + a 



^ , , 1 f mr \^3-2y-^ ^ /TOr + a. /2, 

~3j2(mr)mr— W ■ — [ W ] W - 

m v V + a J r r \ mr V 
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In terms of dimciisionless integration variables 

.8A xii2(0) /2 
A-2 = -I— 5 — 

Aia) ^ 



Aia) 



dx cxp{~x){{x + l)(-j2{x) +j'2{x)x) 



xx{ 



x 



X + a 



X 



X + a 



12 



X + a 



X + a 



X 



-1) 



X + a 



X 



)]} 



-9j2(a:;)i 



X 



X + a 



[1- 



Expanding we find that 



A{0) 
A\0) 



dx exp{—x){x + l){—j2{x) + j'2{x)x)x 
dxexp{-x){{x + l)(-j2(x) +j'2{x)x) 



1 



12 



1 , 1 



+ ^(^)]}-9,2(x)(l)h^]) 

d.Texp(-a;)((a; + l)(-j2(2:) + j2(a;)x) 



,1 6 , 



9j2(a;) 



(F43) 



2 ' 5a;2-' 2x 

Since ^'(0) is finite we can neglect its contribution to the lowest order in a. Fory^(O) we use results below of 



and thus 



.Stt^Q! d 
Ara^ dr 



(F44) 



(F45) 



(F46) 



Similarly, 
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The rate given by Eq. (3.104) is 
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320 

in which we have assumed that 



we have therefore 



r(3Po ^ 27) _ 15 



r(3P2 ^27) 4 ■ 

Our amplitudes and rates for the •^Pq and decays agree with those first computed in [39]. 
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a. Summation Evaluations 



Using Eq.(F17) and appropriate powers of the dummy variable z we obtain 
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For I = j ~ 1 we require the amplitudes 

3™^ dr^ 
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Using once again Eq.(F17) and from above the result K2.1 = 1 as well as 
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and so 



For the next sum we need 
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And finally we compute 
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and 



Now, in addition to the result 
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The final sum to evaluate is 



and so 



Summarizing, 
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